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Abstract 

The intrinsic 4-point coupling, defined in terms of a truncated 4-point func- 
tion at zero momentum, provides a well-establislied measure for the interac- 
tion strength of a QFT. We show that this couphng can be computed non- 
perturbatively and to high accuracy from the form factors of an (integrable) 
QFT. The technique is illustrated and tested with the Ising model, the XY- 
model and the 0(3) nonlinear sigma- model. The results are compared to those 
from high precision lattice simulations. 



1. Introduction 



The intrinsic coupling ^f^, also sometimes called 'physical' or 'renormalized' coupling, is a 
quantity of great interest in a Quantum Field Theory (QFT), especially for scalar fields. 
In some cases, such as the theories, its vanishing implies actually that the theory is 
trivial in the sense that the higher correlation functions of the scalar field can be written 
as sums of products of two point functions, as in a free theory [^. On the other hand, a 
non- vanishing (yf^ is not sufficient to assure the non-triviality of a theory; it only assures 
that a certain four point vertex function does not vanish identically, but does not exclude 
that it vanishes on shell. 

Aside from that, g^i is certainly a renormalization group invariant and a characteristic 
physical quantity of a field theory. In particular it can be used to check the equivalence 
or nonequivalence of different definitions of theories; this will be our main theme in this 
paper, is proportional to the connected - sometimes called truncated - four point 
function at zero momentum, divided by the square of the zero momentum two point 
function and appropriate powers of the mass gap to make it dimensionless; details will be 
given in the body of the paper. 

We are dealing in this article with two main approaches to the construction of a QFT. 
The first one starts from a suitably regularized functional integral and then removes 
the regularization in a controlled way. This is a rather general procedure usable for a 
wide variety of models; it has been successfully employed to construct QFTs in 2 and 
3 dimensions obeying all the required axioms (see for instance 0). Here we will make 
use of a Euclidean space-time lattice gulator. Removal of the regularization, i.e. 

taking the continuum limit in a lattice theory requires the existence of a second order 
phase transition point at which the characteristic length (correlation length) of the model 
diverges. This approach raises the problem of 'universality', i.e. the question whether 
different regularizations yield the same QFT after the regulator has been removed. 

The other approach studied here is applicable to a large class of so-called integrable 
models. It is not based on a Lagrangian, rather the dynamics is specified in terms of a 
postulated exact 'bootstrap' S-matrix, supposed to enjoy a factorization property that 
allows to express all S-matrix elements in terms of the two-particle S-matrix [Q. In 
physical terms this property is linked to the existence of an infinite number of conservation 
laws and the absence of particle production. The postulated S-matrices are then used to 
set up a system of recursive functional equations for the form factors; solving this system 
one can in principle compute exactly all the form factors, in other words continue the 
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S-matrix off the mass shell Once the form factors are known, one can express the 

correlation functions of the basic fields as well as other (composite) operators by inserting 
complete sets of scattering states between them. This gives the correlation functions as - 
hopefully rapidly converging - infinite series of convolution products of form factors. In 
particular in this way one can express the intrinsic coupling in terms of the form factors. 

In both approaches, in principle one has to verify in the end that the axioms of a QFT 
hold. In the lattice approach with a reflection positive action, such as the standard nearest 
neighbor action, essentially the only nontrivial question besides the existence of a critical 
point concerns the restoration of Euclidean (Poincare) invariance in the continuum limit. 
In the form factor approach it is less obvious whether the axioms hold, in particular for 
the form factor expansion of multi-point correlation functions. There exists, however, a 
formal proof (disregarding convergence aspects) of locality and it is hoped, of course, 
that the other field theoretic axioms will hold as well, because the construction is to a 
large extent inspired by them. It is also not clear from first principles - though in practice 
there are very natural guesses -, which field in one construction should be identified with 
which form factor sequence in the other. In any case, even assuming that all the axioms 
hold in both constructions and that one has correctly identified the fields, it is a nontrivial 
question whether the two approaches define the same theory, and in particular whether 
they give the same value for g^. 

In the lattice approach the intrinsic coupling of the two-dimensional 0(n) models we are 
discussing here has been widely studied, both by Monte Carlo simulations 0, H, and by 
various expansions in a small parameter [10, 11, |1^, 13]. For a more precise comparison 
with the form factor approach, we also carried out our own high precision Monte Carlo 
simulations which are reported in this paper. 

In the form factor approach the series for g^^ being a low energy quantity, is expected to 
converge very rapidly. Our results give every indication that these hopes are fully justified, 
though the actual computations turn out to be surprisingly intricate. In the present study 
we want to develop this computational framework, outline the computations and compare 
their results, where possible, to those obtained numerically from the lattice approach by 
the different methods mentioned above. 

Remarkably in all the examples considered the first non-trivial term in the series, which 
contains only one and two particle intermediate states, appears to give about 98% of the 
full answer (!). Moreover for this dominant contribution a general model-independent 
expression in terms of the 1-and 3-particle form factors and the derivative of the S-matrix 
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can be obtained. 

In this paper we discuss three models which can be viewed as the 0{n) nonhnear sigma- 
models for n = 1,2,3. Though formally members of the 0(n) series of nonlinear sigma- 
models, the physics of these systems, their form factor description, and not the least our 
motivation to study them is very different: The n = 1 case is just the massive continuum 
limit of the Ising model. Here the spin form factors are very simple and we were able to 
push the computation of the series up to all terms with a total particle number (summed 
over the three intermediate states) of less or equal 8. The extremely rapid decay of the 
terms is manifest and we use the observed pattern as a guideline for the other systems. 
The final result amounts to a determination of (7r with an estimated precision of better 
than 0.001%. 

The n = 2 case is better known as the XY-model. Here we rely on a bootstrap description 



of the model, to which we hope to return in more detail elsewhere \T^]. Not all the form 
factors are known explicitly, but the specific version of the 3-particle spin form factor 
needed for the dominant contribution can be found by elementary techniques. We compare 
this leading order result with that obtained by lattice techniques and find reasonable 
agreement, which can be taken as support for the proposed bootstrap description. 

Finally the n = 3 model is the first with a nonabelian symmetry group. The evaluation 
of here is in part motivated by the controversy about the absence or presence of a 
Kosterlitz-Thouless type phase transition; see |^ for a more thorough discussion. 

Let us remark that the form factor bootstrap has also been applied to the computation of 
in the sinh-Gordon model; in this model the intrinsic coupling is especially interesting 
because of its relevance to the issue of "triviality" versus "weak-strong-duality". For 



details see the accompanying paper |15 



The article is organized as follows. In the next section we describe the form factor con- 
struction of Green's functions in terms of form factors generally and derive the formula for 
the dominant contribution to the coupling. Further we prepare the ground for the com- 
putation of the sub-leading terms in the specific models. We then give a few generalities 
about the Monte Carlo simulations, and move on to discuss the three models as outlined 
above one by one in more detail, comparing the results of the form factor construction to 
those obtained by the the lattice definition of the models; for the latter the values of g^^ 
are estimated by high temperature expansions as well as Monte Carlo simulations. 
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2. Construction of Green functions in terms of form factors 

In this section we will consider a general massive QFT described in terms of its generalized 
form factor sequences by which we mean matrix elements of local operators between 
physical states. We will restrict our attention to the case ofd = 2 dimensions (although the 
extension to arbitrary d is often straightforward). The application of the representation 
to the integrable models where the form factors are explicitly known will be the subject 
of the next chapter. 

2.1 Generalities 

Our first goal is to construct the Euclidean correlation functions (Schwinger functions) 
from the generalized form factors. The Schwinger functions are convenient because they 
have simpler properties than the Wightman functions and also because it facilitates the 
comparison with lattice results later. For points Xk G IR^, k = 1, . . . ,L, we denote by 
{xki,Xk2) their components and by LXk = {—ixk2,Xki) a Wick rotated version. For defi- 
niteness we will consider here correlation functions of n scalar fields $°(x), a = 1, . . . ,n 
(the generalization to other types of fields is straightforward). Then 

S"^{x,,...,xl) = {^'''{x,)...^''^{xl)) , 

S'''-''^{x,,...,XL)=W'''-''^{LX,,...,iXL) for X,2>...>XL2- (2.1) 

The first equation is the usual operator interpretation of the Schwinger functions. The 
second equation ( p.l|) then indicates the relation of the Schwinger function to the cor- 
responding Wightman function for points (zi, . . . , Zl) = {tXi, . . . , lxl) in the "primitive 
tube" of analyticity.[| Outside the primitive tube the Schwinger functions can in principle 
likewise be obtained from the Wightman functions by analytic continuation and are then 
found to be completely symmetric in all variables. In a form factor expansion however 
the primitive domain is preferred in that only there the convergence of the momentum 
space integrals is manifest through exponential damping factors (c.f. below). We thus 
mimic the effect of the analytic continuation by performing the symmetrization by hand 

^'^-"^ {x,,...,xl)= Yl Se''-"'' i^si, . . . , x,l) , 

S^e'-''^ixu...,XL) := e(a;i, . . . , xl)W''^-''' (cx,, . . . , lxl) , (2.2) 

*We use the signature (+, — ) for (complexified) Minkowski space in which case (zi, . . . , zl) is in the 
primitive tube iff ~lm{zk — Zk+i) G V'^ , A: = f , . . . , L — f , where is the forward hght cone. 
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where Q{xi, . . . ,Xl) is a generalized step function that vanishes unless Xu > • • • > Xl2 
holds and the sum is over all elements of the permutation group S^. The functions 
Sq"'"'^{xi, . . . , xl) are expected to have a convergent expansion in terms of form factors 
in the interior of their support region (as well as for certain points on the boundary). The 
cases of interest here are L — 2 and L — 4. Formally inserting a resolution of the identity 
in terms of asymptotic multi-particle states 1 = J2m Im) (ml obtains 

S^"''{xi,X2) = e(xi,X2) ^e-(^^-^^)^-^2Le^(^i-=^2)iPH,^0|$»i(0)|m)(m|$"^(0)|0) , (2.3) 

m 

and 

S^"'^''"'\xi,X2,xs,X4) = e{xi,X2,xs,x^) ^ e-^^^-^^^^^^e^^^i-^^^i^^ 

k,l,m 

^-{x2-X3)2Ei_^i{x2-X'i)i_Pi_ ^-(x3-X4)2Ern^i{x3-X4)iPrn 

(0|$"i(0)|i)(^|$"2(0)|i) (i|$"3(0)|m)(m|$"^(0)|0) . (2.4) 

The states \rri) are assumed to be improper eigenstatcs of the momentum operator P^, and 
Em, Pm denote the eigenvalues of Pq, Pi on \m), respectively. To write down an explicit 
parameterization of the complete set of states \rn) requires of course the full knowledge 
of the spectrum of stable particles. This is a basic input assumption for the integrable 
models dealt with in the next section. Here for simplicity of notation we will consider the 
case where there is only one multiplet of stable particle states of mass M. An exphcit 
parameterization will then be given in subsection 2.3. 

We introduce their (dimensionless) Fourier transforms V by 

{2nfS^^\h + ... + kL) M-^(^-^) F'^i-«^(A;i, . . . , A;^) 

- J d^xi . . . d'^XL S^'-""^ (xi, . . . , Xl) e^('=i^i+-+*^^^^) , (2.5) 

taking into account the translation invariance of 5*0. The Fourier transform of the full 
Schwinger function is then obtained by symmetrization 

S^^-^^{ku ...,kL) = {27rf 6^^\ki + ... + kL)M-^^'^'^^ . . . , ksL) , (2.6) 

seSL 

whereon rotational invariance gets restored. The desired representation of the two and 
four point functions in terms of form factors is given by 

v^^^-(k,,k2) = ^y^^"^(fci,A;2), 

m 
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where 



V-^^-^^^^{k„h,h,h) = J2y^''''''ik,,h,h,h), (2.7) 

k,l,m 



— Era- iki2 



Ek - iki2 Ei_ - iki2 - ik22 + ik42 
(0|<l>''^(0)|i)(^|<l>"2(0)|i) (i|<l>"^^(0)|m)(m|$«4(0)|0) , (2.8) 

with the understanding that the sum of the momenta kj vanishes. Further we denote by 
V^^^^iki, fcz) and ^fc'l^'"'"''(A;i, fca, ks, k^) the quantities (U) with the integrations over the 
rapidities performed, the measure being inherited from Eq. ( |2.25[ ) below. 

The key assumption of the form factor approach in this context is that the matrix ele- 



ments in ( |2.8| ) can be computed exactly via solutions of a recursive system of functional 



equations, the so-called form factor equations or Smirnov axioms. Symbolically 

il\^-{Q)\rn) ^ ^,';...,^,^...,,„(^i,...,^,|0i,...,^^)=:^^^(^|0). (2.9) 

The rhs, for which we shall often use the indicated shorthand notation, is called a gener- 
alized form factor, the special case with either / = or m = are the form factors proper. 
The form factors are meromorphic functions in the rapidities, while the generalized form 
factors are distributions. The form factors can be computed, at least in principle, as 
solutions of the before mentioned system of functional equations. The generalized form 
factors can then be obtained from them by means of an explicit, though cumbersome, 
combinatorial formula. We shall later just state the special cases of this formula required. 
A discussion of the general formula can e.g. be found in the appendix of ||16 |. 



Implicit in the products of matrix elements in (|2.8| ) of course are appropriate index con- 
tractions. For definiteness let us note them explicitly 

(0|$'^(0)|m)(rn|$^0)|0) ^ C''(0) , 

m\m){m\m^mm\m){rn\^\m) — 4t;Viei^), (2.10) 

where 

A 

4"tviei^)= E n{^)^\-B{^^\onTc{e\o)ri-{e^)- (2.11) 
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Here = (a^, . . . ,ai), u'^ = {ook, ■ ■ ■ ,1^1), etc. The construction is such that I^i^) is 
a completely symmetric function in ^ = (^^i, . . . , 6^)- Similarly Ikimi^l^l^) symmetric 
in each of the sets of variables uj = {lji, . . . , a;*;), ^ = ({1, ■ ■ ■ and 9 = (6*1, ... , 9m), 
individually. 

2.2 The intrinsic coupling 

As surveyed in the introduction the intrinsic coupling is defined in terms of the zero 
momentum limit of a connected 4-point function. We may assume that the Schwinger 
functions of the scalar fields with an odd number of arguments vanish; then the connected 
L — 2,4 Schwinger functions of interest here are 



- >5"i"3(A;i,A;3)^"^"^(A;2,A;4) --S"i"^(A;i,A;4)-S"^"^(A;2,A;3)(2.12) 



Making explicit the overall delta-functions arising from translational invariance we intro- 
duce the Green functions by 



where the constraint ki -\- . . . -\- ki — Q in the arguments of Ci - ^i will always be under- 
stood. 

In the following we will now assume that the theory is 0{n) invariant and thus for the 
2-point function we can write 



(2.13) 



G''"'^{k,-k)^S''"''G{k). 



(2.14) 



The intrinsic coupling is then defined by 



gn = -^f■ 



M2 1 

G(0j2^ 



-^G-^^(0, 0,0,0), 



0,6 



(2.15) 



where we leave the choice of positive constant J\f for later. 
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Performing the symmetrization ( ^.21 ) and the Fourier transform one recovers the famihar 
expression for G{k) in terms of the spectral density 

G{k)= / dfipifi)^—-, (2.16) 

Jo jjj -]- K 

where 

p(/i) = E ^(/^ - - P^) 4vri5^5(Pj- E(0|'^"(0)Im)H*"(0)|0) . (2.17) 



In order to compute S'^'^'^^'^'^""^{ki, k2, k^, k^) the symmetrized sum ( p.6|) , ( p7|) has to be 
performed. For reasons that will become clear immediately we first single out the partial 
sum with / = 0. Taking into account the ^4 permutations one finds 

{2nf5^^\ki + k2 + k, + k,)M-^ J2 E VkoT^^'ik,,, k,2, ks3, k^,) 

k,m s£<S4 



S^'^'iku k2)S''''''ik3, ki) + ~S^^^^{ki, ks)S'''''^ik2, ki) + ~S^'^\ki, A;4)^'^^"^(fc2, k. 



3 J 



+ {2'nf5^^\ki + k2 + k3 + ki)M-^ J2 ^''''"'"''"'''iksi, A;,2, ^,3, ks,) . (2.18) 



Here 



l^"^"^"^"*(fci, k2, fcg, k^) = -^5iku + k2i)M^ 5''""'5'''''^Hik^, k2)G{k^) , (2.19) 

with 

H(k,,k.)=l (2,20) 

In obtaining ( p.l8|) we defined the second denominator in (|2.8|) for / = with the ie 
prescription as: —i{ki2 + A;22 + ^e)- Here and later the distributional identity 

V tTTO[X) , 



X + ie X 

will be heavily used, where V is the Principal Value prescription. 

One observes that the first three terms in ( p.l8|) are precisely the ones removed by the 
definition of the connected 4-point function. Remarkably there is a remainder, the fl term, 
which is present even in the free theory. Typically the spectral densities are decreasing or 
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bounded by a constant as — >■ oo. The functions G{k^ and H{ki, ^2) are then regular at 
ki = 0. Inserting finally ( p.l8| ) into ( p.6| ), ( |2.7] ) one obtains for the Green function ( p.l3| ) 

= E E ^^""''^'""(^si, A;.4) + E n'''"''"''"'''{ksi, h2, ks3, K,) . (2.21) 

On general grounds one expects the vertex function to be real analytic. In particular 
there must also be terms involving the delta function in the Vkim above which cancel 
those of the fi-term (we will demonstrate this explicitly in the computation of V121 in 
appendix A). Thus, provided the coupling is well-defined (finite) at all, the result will be 
independent of the way the zero momentum limit is taken. It is therefore desirable to 
find a convenient limiting procedure that simplifies the computation. To this end we first 
observe that in (|2.8| ) the kji and kj2 components enter asymmetrically. In particular as 
long as the intermediate state is not the vacuum (i.e. / 7^ in the second formula, and 
recalling that we are assuming that none of the operators involved has a non-zero vacuum 
expectation value) one can put kj2 = 0, j = 1, 2, 3, 4. We now compute the 4-point vertex 
function at zero momentum through the limiting procedure: 



G'"!"^"^"^ (0,0, 0,0) = lim G''^''^''^''^{ki,k2,k3,ki] 



kj = {MshKj,0) 



(2.22) 



where J2'j=i sh.Kj = 0, and the limit is taken such that 



7^ l^jl for i j , and — Kj\ 7^ — for distinct pairs. 



(2.23) 



In view of ( p.l9|) it is clear that the limit prescription in (|2.22|) , which we will use in the 
following, has just been designed such that the ^2 term does not have to be considered in 
computation of the coupling. 

Before embarking on further computations let us comment on a few structural issues. On 
physical grounds one expects the intrinsic coupling to be both finite (in a theory with a 
mass gap) and positive (for Af > 0) when the interaction is repulsive. Mathematically 
however it is a quite challenging problem to actually prove this, whatever non-perturbative 
definition of the theory one adopts. In the context of constructive (lattice) QFT such re- 
sults seem to be available only for a single phase $2 theory (see e.g for a survey). 
In the present context we wish to define the theory strictly terms of its form factors. 
Mathematically speaking one should then try to prove in particular that the right hand 
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side of ( p^.2lD defines a real analytic function. For the dominant low particle contribu- 
tions we demonstrate in appendix A explicitly that all non-analytic (e.g. distributional) 
terms indeed cancel out. We have not attempted to prove this in general, nor can we 
estimate the rate of convergence of the sums in ( 2.21 ) on general grounds. In all the 
examples considered later however the series appears to be rapidly convergent; the terms 
are alternating in sign and decrease in magnitude very quickly with increasing particle 
numbers. 



2.3 State parameterization 

Here we assume that the single particle spectrum consists only of an 0(n) vector multiplet 
of mass M. The one particle states |a, a) are thus specified by an internal "isospin" label 
a and the rapidity a (i.e. the spatial momentum of the state is p = Msha). The states 
are normalized according to 



(a, /?) = 47r(5afe5(a — {3) . 



(2.24) 



The condensed notation for the sum over states now becomes 



|l?l)(l?2;| 



- |o)(o| + E E 

m=l ai,...a 

°° d^i r^i d^2 /■''—I d^, 

-oo 



m 

oo AtC 



ai,6i; . . .; am, OmY'^ "(c^i, ^i; • • • ; Om, Om\ ■ (2.25) 



It is often convenient (for a fixed m) to perform the change of variables 



«i = - ^i+i 5 J = 1, • • • , - 1 



A = - In 



(2.26) 



since in terms of these variables the total energy and momentum of the states take on a 
simpler form: 



{Em, Pm) 



M ch^j , M sh^j' ) = (m(™) {u) chA, M^"^) {u) shA 



(2.27) 



where the eigenvalues Mm = jE^i — Pm of ths mass operator are given by 



M^"") (u) = M m + 2 ^ ch{ui + ... + uj^i) 

i<j 



1/2 



(2.28) 
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Correspondingly the integration measures in (p^.25D above are replaced by 



°o d^i fSi d^2 f^"^-^ d^m d^-^u dA 



oo 



(2 29) 

in 7-00 in "' J-oa 4:71 Jo (47r)'^-i J-oc An ' 

For later reference we also display the inverse transformation 

Oj = Uj + . . . + Um-l + Mm + A , j = 1, . . . , m , 

where := In , , '^-i , • (2-30) 



m — 1 



2.4 The two point function 



The spectral function (|2.171 ) appearing in the representation of the two point function can 



be written as a sum of contributions of fixed particle number m 

P(/i)= E P^"^-"), (2.31) 

0<m odd 

where only odd numbers of particles contribute due to our assumption that the fields $° 
are parity odd. We normalize the fields by 

(0|$'^(0)|6,a) = 5,% (2.32) 
rendering the 1-particle contribution to the spectral density simply 

p(i)(^) = 5(^-M). (2.33) 
The m > 3-particle contribution to the spin spectral function (|2.17| ) is given by 

p(™)(/i)=y^ J^-^6{^^-M("^\u))Uu), (2.34) 

with 

Im{n) ■■=-!: E \^L.aJ0u---,0m)\', (2.35) 

a ai,...,am 
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which equals I}^{9) under the integral. The function featuring here corresponds to 
the matrix element of $° between vacuum and an m-particle in-state as in ( p.9| ) 



K...aJ^^^ • • • , ^m) = (0|$"(0)|ai, e^-...- a„, , 9^<...<9^. (2.36) 

The inverse 2-point function has a low momentum expansion of the form 

G{k)-^ = Z^ '[Ml + fc2 + 0(A;")] , (2.37) 



with 



2 

Ml = M^^ , = ^ , (2.38) 

02 02 



where 72, S2 are spectral moments: 

72 = r %p{f^) , 62 = M' r ^p(/x) . (2.39) 



2.5 The intrinsic coupling revisited 



In ( p.l5|) we left open the choice of the normalization constant M because for different 
models different choices are convenient. In analytical and numerical lattice computations 
(at fixed cutoff) it is often easier to compute the second moment mass Mr instead of 
the (exponential) spectral mass M (in lattice units). For ease of comparison with these 
techniques we thus choose M = M'l/M'^, i.e. we define the intrinsic coupling by 

/lf2 1 

3- = -7^2- E ^'"'^"'(0, 0, 0, 0) . (2.40) 
Using 0(n) symmetry it follows 

Qaia2asa4(^Q^ 0, 0, 0) = M-^-f^iS"'"^ 5"'"'^ + 5'^i'^3^a2a4 ^ §a,ai^a2as^ ^ ^2.41) 

and hence we can write ( 2.401 ) as 

71 + 2 74 , , 

9n = X • 2.42 

n 72^2 
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These spectral moments have, corresponding to the decomposition ( |2.31D , an expansion 
in contributions arising from states with a fixed (odd) number of particles 



72 = 1+ 72;m, ^2 = 1+ E ■ (2.43) 

3>m odd 3>m odd 



Similarly, corresponding to the sum in (|2.21|) we have 



74 = E 74;fc«m , 74;fc/m = 74;mZfc , (2.44) 

k,l>0,m 

where the sum goes over odd integers k, m and positive even integers /. To avoid writing 
many 0{n) indices we will use 



lA;klm = ^Lim^^^O E ^klm{tisl, K.s2, t^sS, ^^m) , (2.45) 



where 

1111, 



Vklm{f^l, f^2, %5 f^i) — Vklm (^l5 ^2, ^3? ^4 



(2.46) 

kj={MshKj,0) 



and the symbol Lim above means taking the limit kj — > with the kj satisfying J2j shuj 
and the constraints in Eq. ( p.23| ). 
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3. The nonlinear 0(n) sigma- models 



As outlined in the introduction the form factor bootstrap (FFB) construction of an in- 
tegrable quantum field theory starts from postulates of the on shell properties of the 
theory. By integrable here it is meant that the theory has an infinite set of conserved 
charges which entail that there is no particle production. This property usually is char- 
acteristic of non-relativistic quantum mechanics, remarkably here it holds for relativistic 
quantum field theories (QFTs) (assuming that the FFB approach does indeed define a 
QFT). In 4-dimensions absence of particle production implies that the theory is free but 
in two dimensions this is not so. In addition to the absence of particle production, one 
postulates the spectrum of stable particle states and their 2-particle S-matrix which has 
to satisfy the so-called Yang-Baxter (or factorization) equation ( |A.5| ). 

In principle one could proceed without reference to a Lagrangian, but often contact to a 
Lagrangian description is desirable. Thus typically postulates of specific S-matrices are 
motivated by studies of associated Lagrangian QFTs. Unfortunately in most cases one 
cannot solve the QFTs to the extent necessary to really derive the candidate S-matrix, 
rather one has patches of partial information. This is in particular the case for the 0(n) 
nonlinear sigma models formally described by a set of spin fields a", a = 1 . . .n > 2, with 
the constraint = 1 and Lagrangian density oc (9^cr)^. There is a wealth of information 
on these models which will be recalled when we study the various cases in the following 
sections, and for an overview we refer the reader to our previous paper 0. In particular 
the spectrum of stable particles is thought to consist of an 0(n) vector multiplet of mass 
M without further bound states (i.e of the form of the spectrum considered in subsection 
2.3). The S-matrix element (for n > 2) has the decomposition 

Sab;cd{d) = SabScd + (^2(6) Sac^bd + O^iO) 5adhc , (3.1) 

where the center of mass energy is given by -y/i = 2Mch6'/2. 

Classically the theories have an infinite set of local and non-local conserved charges. One 
can argue that there are no anomalies which obstruct the existence of such charges in 
the quantum theory. In the case of the non-local charges for n > 3 the construction 
of Liischer |jT^ is closely connected to the usual perturbative renormalizability and the 



(perturbative) asymptotic freedom of the model. Knowledge of the action of the non-local 
charges on the asymptotic states then restricts the S-matrix to the form postulated by 
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Zamolodchikov and Zamolodchikov for n > 3 

-2Tiie S2{9) 



am 



[i-K -9) (n - 2)9 - 2m 



= -2ni 



S2{9) 



(n - 2)9 - 2m ' 

i.e. the invariant amplitudes are all given in terms of one amplitude which we have chosen 
here to be the invariant amplitude S2{9) in the symmetric traceless ("isospin 2") channel. 
The amplitude ^2(6') is off-hand determined only up to so called CDD factors, which were 
initially fixed by selecting the solution with the minimal number of poles and zeros in 
the physical strip. This solution for S2{9) is given by 



= -exp|2i / — sMuj) KJuj)] (3.3) 

with 



^nM= ^l,-.. • (3-4) 

The proposed identification of ( ^.2[ ) - ( p. 41 ) with the S-matrix of the 0(n) sigma- model 
passes several non trivial tests. First, the leading terms of its large n-expansion coincide 
with those obtained in leading orders of a field theoretical large n computation. Second, 
in the determination of the exact M/ A ratio a consistency condition arises when matching 
the results of a perturbative computation against that obtained via the thermodynamic 
Bethe ansatz 14^]. This consistency condition is also sensitive to the CDD factor; the 



minimal bootstrap solution ( ^.2[ ) - ( |3.4| ) passes the test. 

We note that the above formulae have a smooth n ^ 2 limit. A study of the possible 
relation of the so defined FFB 0(2) model to the continuum limit of the lattice XY model 
(from the massive phase) will be the topic of a future publication Jl^ . 

Further we remark that the S-matrix for the case n = 1 (Ising model) can also be written 
in the form ( p.l| ) by setting 

ai{9) = <y2{9) = Q , a;{9) = -I , n=l. (3.5) 

The representation ( p.l| ) is of course redundant in this case, but it does allow us in 
the following to discuss all n > 1 simultaneously. For example in all cases we have an 
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expansion at low energies of the form 

Sab;cd{0) = -6adSbc + ^0Dab■,cd + 0(6^) , (3.6) 

in sharp contrast to a weak perturbation of a free field theory. 

Having all the on-shell information covers all the physical information on the theory one 
observes from scattering of the stable particles, but off shell information is being explored 
if the system is probed by external sources weakly coupled to local operators with given 
quantum numbers. 



3.1 Derivation of the leading term in the FF expansion of 

For the leading 1-2-1 particle contribution to a general model-independent expression 
can be given in terms of the derivative the S-matrix and the 3-particle form factor. For 
notational reasons we restrict attention here to the 0{n) models considered later in more 
detail. The extension to a general integrable QFT without bound states is described in 
appendix A. For the 0(n) models the formula reads 



3 



74;121 = 4i X! 




Mu)M-u)-^]. (3.7) 

ch M 



Here /c(^) is a particular version of the 3-particle form factor T^^^{6i,62,6'i) of the lo- 
cal field supposed to correspond to the renormalized spin field o"^ in a Lagrangian 
construction. Explicitly 

W):=J'l,{m,e,-e). (3.8) 

In order to derive ( ^.7| ) consider first more generally the (1, Z, 1) contribution in ( |2.21| ) with 
/ > 2. Using (|2.8|) and switching to the explicit notation introduced in subsection 2.3 one 
can perform the integrals over the rapidities of the '1' particles. Then one decomposes the 
rapidity measure for the intermediate 7' particle contribution according to ( |2.26| ). Using 
( |2.29| ) the A integration can be performed and by means of (|2.3CI| ) one arrives at 



vM-un.,n,,n,)= 2chVchV/!i (i^^E^^^ '^"^""^'''"^^ 

1 r d'-^M 1 m2 



f d u 1 M 



8ch\ich2K4 7 (47r)'-ich^A, M(0(m 
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where Im{—K,i\6\KA) := IIII^{—Ki\6\k,4) is a product of generalized form factors as in 
( |2.1CI| ), ( p.ll| ). Explicitly the correspondence to the matrix elements is 



bi,...,bi 

x'-{bi,ei;...;bi,ei\<!>\0)\l,K^). (3.10) 
In the first expression we introduced the notation 

S{9„...,9l) = 5(j2sh9^ , (3.11) 
in the second one A* is defined by 

M 

shA* = - j^^i)^^^ (sh/ti + shK2) . (3.12) 

As remarked before the generalized form factors can be expressed in terms of form factors 
of the same operator and delta distributions by an explicit combinatorial formula. We 
shall usually just display the specific version needed. A discussion of the general formula 
can be found in the appendix of |T^. For the generalized form factor entering ( p. 10 ) the 



formula reads 

{ai,-K,\^-'{0)\h,9,;b2,92y 

b2;a2ai (^1 92)6{k, + 92). (3.13) 



= J'T,,,,2i-^l+^7r-^e,9^,92) 

1>6»2 



Substituting this in Eq. ( |3.9| ) we obtain 

1 r , 1 



647rch^Kich^ cha2(chai + cha2) 

xj^lj.y{in — Ki — ie, ai, a2)J^lxy{iTr — k^ — ie, —ai, —02) 
Stt 



chK3(chK3 + chK4) 



J^lii{iTT — Ki — ie, K4, K3) (3.14) 



+ Ch.2(ch!"+Ch.2)-^'"^'" " ~ "^^^ • 

Here we used the simplifications discussed above and the real analyticity property ( |3.16| e) 
below. Moreover, we changed the integration variable from the difference of the two 
rapidities to one of the rapidities (ai). The other rapidity {d2) is then the solution of the 
transcendental equation 

shai + shQ;2 + shui + shK2 = (3.15) 
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and is an analytic function of ai. There are no contributions from terms involving delta- 
functions such as S{k,i + K4) appearing in J121 since we are taking the limit ( |2.22| ) where 
these delta-functions vanish. (These terms are however crucial to cancel corresponding 
singularities in the Q term; c.f. appendix A). 

The form factors appearing in (|3.14j) obey a system of functional equations which allow 



one to further simplify the expression. Let us recall these equations in the form relevant 
to the three-particle form factor jF^^^(a, /5, 7) in the 0(?7,) model. 

P, 7) = 5be;,x.(/5 - 7)^a.,(«, 7, P) , (3.16a) 

P, 7) = -^1(7 + 2m, a, (3) , (3.16b) 

/?, 7) = ^1(« + A, /? + A, 7 + A) , (3. 16c) 

^l(a, A 7) = ^l(-7, -a) , (3.16d) 

7)]* = , -7*) • (3.16e) 



Here the S-matrix appearing in the exchange axiom ( |3.16| a) is the 0(n) S-matrix ( p.l| ). 



( |3.16| d) and ( p.l6| e) express the parity invariance and real analyticity property of the 



form factors, respectively. The homogeneous axioms ( |3.16D are supplemented by the 
inhomogeneous residue equation 



lim {a- [3- in)J^^^^{a, P, 7) = 2i{5ab5cd - Sbc;ad{P - 7)} • (3-17) 



We now take advantage of the analytic properties of the form factors and change the 
integration contour in (|3.14|) from the real axis to a curve C which is arbitrary except that 
it has to stay within the 'physical strip' < Imai < 7r/2. Along this contour we can put 
e = and also the limit can safely be taken. The integrated part of ( |3.14|) then 

simplifies to 



1287r Jc ch^a 
where we introduced the shorthands 



viS = 7^ f 4^Ma)M-a) , (3.18) 



/i\c(") =: ^bcfb{a) (no sum) . (3.19) 
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Of course, one has to take into account the contribution of those singular points of the 
integrand that get crossed when deforming the contour of integration. There are two such 
singular points: 



ai = K4 + ie and 0^2 = —^1 — 'i^ , (3.20) 



which never coincide if ( ^.23 ) holds. 



Applying Cauchy's theorem one can evaluate the contribution from the first singular point 
using the residue axiom ( p.l7|) . This gives 



16ch Kich K4chQ;2(chai + chQ;2 



{•^iii(^^ — Ki — ie, «!, 0:2) — -^iii(^7r — fti — ie, 02, a;i)| 



where in the second term we also used the exchange axiom ( |3.16| a). After taking the 
limit e —>■ 0, which is possible if ( |2.23 ) holds, the contribution of the first singular point 
becomes 

-—2 -2 — — r-^{-^ni(^^ - ^3, /^4) - - Ki, Ki, K3)} . 

loch Kich K^chn^yaiK^ + chK4j 

The contribution of the second singularity is similar: 
1 



[j^llSir - K4, K2, Ki) - J^i\i(ivr - Ki, Ki, K2)}. 



16ch^fi:ich^K4chK2(chKi + chK2) 

Putting together the contribution of the singular points and the last two terms of (|3.14|) 
the non-integrated contribution can be written as 



8ch^Kich^K4chK3(chK3 + chK4) 



{j^uiii-rr - Ki, K^, K4) + J^ni(ivr - Ki, k^, K3)} 



where = indicates equality after the symmetrization over the elements of the permutation 
group 1S4 has been carried out. 

We now use the Smirnov axioms (|3.16|) and ( p.l7|) to simplify the non-integrated part in 



the (symmetrized) — *■ limit. It is convenient to first introduce the reduced form factor 
Oabci^,P,l) by 

J-l(a, (3, 7) = T,{a, (3, 7)^1(«, A 7) • (3.21) 
Here and in the following we set 

TNi^^,...,^N)■.= n (3-22) 

l<j<j<Af 
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where T is basically the tanh-function T{6) := tanh^/2. Note T{6) has a simple pole 
at ^ = iir, T{m - B) = -T{iTi + 9) = -2/6 + O(^), and a simple zero at ^ = 0. The 
advantage of the representation (|3.21|) is that the singularities are carried by the tanh 
factors and the reduced form factor ^^^^ is analytic everywhere in the physical strip. In 
particular, for small a,/? and 7 it can be expanded as 



^KzTT + a, /3, 7) = JL + (« - 7Xc + (/5 - 7)^L + 



(3.23) 



where the dots stand for terms higher order in a, (3 and 7. We can compute the constant 
tensors appearing in the expansion ( p. 23 ) using the form factors equations. From the 
residue axiom ( |3.17|) we can immediately fix 



Jibe = i{^ab^cd + ^achd) 



^abc + ^abc ~ ^hc-ad ■ 



(3.24) 



To determine the expansion coefficients individually we employ the exchange relation 
(|3.16| a) and find 



^abc — Dhc;ad — Dfyc-da , -^^abc — D^c-da ■ (3.25) 

Using the expansion ( |3.23|) , for small n the non-integrated contribution becomes 



(/) . 1 



K3 — Ki 



4 (ki + K3)(fti + 



(3.26) 



where 



D = D 



11:11 



_^^da,(e) 



d^ 



This can be simplified by noting that upon averaging over the permutations 



(3.27) 



(3.28) 



and similarly 



K4 



1^3 



1^3 



_ 1 

+ K3 2 



Kl + K4 Ki + K3 Ki + /t4 

After this simplification we have for the non-integrated contribution 



(3.29) 



(3.30) 
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and hence the non-integrated part of the leading contribution to the four-point couphng 
eventually becomes 

7^21 = -4/^, (3.31) 

For the Ising model the S-matrix is constant and therefore 74^121 = for n = 1. For n > 2 
we use ( p.3| ) and find 

Ti'm = -- + 8 r dcukr^iuj) . (3.32) 
vr JO 



The integrated part (|3.18| ) (which is in this form rather useful for numerical evaluation) 
can be written in an alternative form using the the residue axiom which implies 



4 , 4? 

Jl + 0{1) = -j{5a,5,a + 5aAd) + 0(1) 



(3.33) 



Using this we can explicitly subtract the singular part in ( 3.18 ) and shift the contour back 
to the real axis. Noting also that the integrand is an even function of a we arrive at 



74:121 



1 \ 1 

— / du< — n— fbiu) fbi—u] 
Svrio Ich^'^^ ^ ' "'2 
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(3.34) 



The extension of the formula (|3.71 ) to general integrable models without bound states is 
described in appendix A. 



3.2 The three particle form factor 



Only the special three-particle form factor fahSf^) iii ( P-19| ) is necessary to compute the 



leading contribution (|3.7| ) to g^. It turns out to obey an autonomous system of functional 
equations (in a single variable) that derives from the form factor equations satisfied by 
.7-'^j^(a, /3, 7). Solving it allows one to compute fh{&) - and hence to evaluate ( |3T7| ) - in 
situations where the general form factors are not known. 

We begin by noting that the functions fthSP) real analytic, i.e. [/a{,c(^)]* = fabci~^*)y 
in the physical strip < lm6 < tt, with simple poles at 6* = and 6' = y. Moreover, 
using (|3.16| b,d) one can easily deduce that it is symmetric in its last two indices. 



fLiO) = fUO) ■ (3.35) 



Using ( |3.16| a) one obtains 



fL{0) = S,c;y.mft(-e), (3.36) 
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and finally combining ( p.l6| a-d) results in 



/aK^vr -e) = s^,y.M SyUm Si.,uo) fL-A^^ + o) 



(3.37) 



These are the consequences of the homogeneous form factor axioms; they are supple- 
mented by the residue equations 



bd 



Res/l(f) 



In view of Eq. ( p. 351 ) we can parameterize / as 



fabM = m5ad5bc + m 5ac5bd + 5abS, 



Jed 



Then the contribution 74^/21 is given by 



(11) 

74:121 



57r JO 



du 



ch^M 



{n - 2)9 + 2m 
{i-n - e){m - 2e)[{n - 2)6 - m][{n - 2)9 - 2m]2 ' 



(3.38) 
(3.39) 



(3.40) 



■ nk{u)k{-u) + 2k{u)l{-u) + 2k{-u)l{u) + Al{u)l{-u) 64 ' 



(3.41) 



(3.42) 



In terms of the two functions k{0) and l{0) Eq. ( p.36| ) can be written as 
k{e) = [s2{2e)+nai{2e)]k{-e)+2ai{2e)l{-e), 

m = sA2e)i{-e), 

while ( p.37|) becomes 

kin,- 6) = [AiAe)k{:LTi + e) + Ai2{e)i{ni + e)]a{e)s2{e)^s2{2e), 

liiTi-e) = [A2i{e)k{m + e) + A22{e)i{n: + e)]a(e)s2{e)^s2{2e). (3.43) 

Here 

a{e) 

and 



(3.44) 



All (0) = {0- in) [2{n - 2)^9^ + {n-2){n- 4)9^71 + {n + 2)97,^ - 2i7c^ 

Ai2i9) = -A{n-2)in9{9-i7r){9 + in), 

A2i{9) = -2{n - 4)t7r^9 , 

A22{0) = Au{-9). (3.45) 
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The matrix A{6) satisfies 

A{e)-' = A{-e)a{e)a{-e) , detA{e) = \ . (3.46) 

The functional equations ( |3.42| ), ( |3.43| ) still contain the transcendental function S2{0). It 
can be eliminated by the following standard procedure. We introduce the function u{6) 
as the unique solution of 

u{9) = S2{9)u{-9), (3.47) 

u{i7i-9) = ~u{i7r + 9), (3.48) 

subject to the normalization condition 

u{i7i -9) = l + 0{9) . (3.49) 
9 

Using the results of Appendix D one can immediately write down the solution 

^(0) = _lT(e)e^W, (3.50) 

where in ( p.3|) of course the kernel Kn^u) defined in ( p^) has to be used. Introducing 

2? 1 

Y{9) = —uiin -9)uiiTT + 9)ui29) , uq = u\0) = — e^^^K (3.51) 
uo 4 

we parameterize /c, / as 

k{9) = Y{9) K{9) and l{9) = Y{9) L{9) . (3.52) 

Rewriting then the functional equations ( |3.42| ), ( |3.43| ) in terms of K and L, the new 
system involves only rational coefficient functions. Explicitly they read 

= U ^aT-^- 7^A[{n-2)9 + iix]{i7: + 29)K{-9)+AiM{-9)], 

[\n — 2)9 — ■i7rj(27r — 29) 

L{9) = L{-9), (3.53) 
and 

K{i'K-9) = [An{9)K{i7T + 9) + Au{9)L{t7T + 9)]a{9), 

L{iix-9) = [A2i{9)K{in + 9) + A22i9)L{m + 9)]a{9), (3.54) 
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respectively. The first equation of ( |3.53| ) can be used to eliminate L{6) in favor of K{d) 
via 



L{e) = -^1 [ct -{n- 2)e]{ni - 2e)K{9) - [in + (n - 2)^](m + 2e)K{-e)} , (3.55) 

and ( p.55| ) also solves the second equation of ( p.53| ). Inserting ( |3.55| ) into ( |3.54| ) results 
in a single linear functional equation for K{6). The normalization of the solution is fixed 
by the residue equations ( |3.39|) . 

We expect that this procedure can be used to compute fb{0) and hence the leading 
contribution to the coupling for all 0{n) models. For the 0(2) model we demonstrate 
this in section 6. 



3.3 Sub-leading contributions 

In order to achieve higher accuracy and to obtain some clue on the rate of convergence of 
the series (|2.44| ) we will compute some of the sub-leading terms as well. It turns out that 
the 1-2-1 term indeed gives the numerically most important contribution to the coupling. 
But based on the computation of the sub-leading terms the numerical result can also be 
endowed with an intrinsic error estimate. Our results indicate that the next important 
contributions to the coupling are (1,2,3) + (3,2,1) and (1,4,1). Its explicit evaluation 
is deferred to appendices C and B. The difficulty in the evaluation lies in the rapidly 
varying nature of the integrands, which have in the multidimensional phase space many 
zeros and (integrable) singularities. To deal with these we have either decomposed the 
integrand into appropriate parts or avoided the singularities by shifting some contours of 
integration into the complex plane. 



The (1, 4, 1) contribution is the / = 4 case of Eq. ( p.9|) and will be evaluated in appendix 
B. Here we prepare the ground for the evaluation of the (1, 2, 3) + (3, 2, 1) terms. More 
generally let us examine the (1, 2, m) + (m, 2, 1) contribution and to this end return to 
(|2.8| ). Performing the internal rapidity integrations one obtains 



2 



1 /■ d^d^ 5(6,6, /«2) 



ch^Ki ml J (47r)2 ch6 + ch6 



X j (4vr)»^ E-ichg, ^-"^(-'^^l^^'^il^) ■ (3-56) 

Next one spells out /12m := 1 12m by inserting the formula expressing the generalized form 
factors in terms of ordinary form factors. Taking advantage of the S-matrix exchange 
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relations many of terms contribute equally upon integration and one ends up with four 
terms 

-jn ^Sin.Ain) (IV) 



with 



H2biaia2...am(P2 + 1^-: A + «7r_, «!, Q;m)^aia2...a^ («1 ' , (3.58) 

(//)/ N 1 V- V- f -.2 {(^1,^2,1^1,1^2) 

~ 8(4.)™-'(™-l)! £ / ' ^ chA + chA 

•^62a2...a..(/^2 + i^-, ^2, am)^6\a2...a^ (A , (^2-, am)* , (3.59) 



(5(/3i,/ 32,ft:i,ft:2) 
+ ch/32 



•^a3...a^(«3, «m)^b\fea3...a^(A' /52, «3, ^m)* , (3.60) 



5(qi. .... n,„,. -K4) 



•^naia2...a^(-'«2 + «7r_, -Ki + «!, Q;m).^aia2...a^ («1> «m)* , (3.61) 

where 7r_ stands for tt — e. All integrals range from — oo to +00. 

Further details of the computation of the 1-2-3 contribution are given in Appendix C. Note 
that for the numerical evaluation of the k -\-m — Q contributions we need the analytic 
expressions for the 5-particle form factor of the spin operator. Unfortunately these are 
at present only known for n — 1 and n — ?>. For the Ising model all the form factors are 
explicitly known and in this case we have also computed the k + l + m — S contributions. 
After a preparatory next section where we discuss the definition and measurement of the 
intrinsic coupling in the lattice regularization, we will discuss the cases n— 1, 2, 3 in turn. 
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4. Lattice computations of g 



In the subsequent sections we will compare the results of the form factor bootstrap cou- 
pling with those obtained from the lattice theory. As noted earlier, in the framework 
of the lattice regularization there are two methods to compute g^^ in the 0{n) models: 
high temperature (= strong coupling) expansions and Monte Carlo simulations. Both 
approaches usually take the standard lattice action on a square lattice 

5 = -/3^a(a;)-a(x + /i), (4.1) 

X,fJ. 

as the starting point, where cr(x) ■ a{x) = J2a'^'^{^)(^"'{^) = 1- 
The lattice definition of 5'r(/?) is as in Eq. 



where all quantities are defined analogously to the continuum theory 

G2{k) = ^EEe^'^(^"(^)^"(0)), (4.3) 

"ax 

-[(a'^na;i)c^"'(a;2))(a"'(x3)fT"^(0)) +2 perms]}, (4.4) 
and ^2 is the second moment correlation length 

^2 = 77^' /i2 = -EE^'(^"(^K(0))- (4.5) 

4^2(0) a X 

The coupling from the lattice regularization is defined as the continuum limit 

g,^ = hm g,i/3) , (4.6) 

where /3c is a critical point where the correlation length diverges (in lattice units). 
Butera and Comi [|1^] have produced long high temperature series for 6*2(0), /i2, and G4 in 



the 0(n) model with standard action, and Pelissetto and Vicari [|13] have reanalyzed these 



26 



series to compute estimates for the intrinsic coupling for n < 4. Similar computations 



have been performed previously by Campostrini et al |]TT 



Our Monte Carlo simulations were of course done on a finite lattice, more precisely a 
square lattice of size L (points) in each direction and periodic boundary conditions, both 
with the standard action ( |4.1| ) and the fixed point action of ref. |]T^ . The infinite volume 
lattice coupling Qr^P) is then obtained as the limit 

gM= \imgR{P,L), (4.7) 

of a finite volume coupling g^i^P, L) which is proportional to Binder's cumulant ul'- 



2 (i^r) 

where E"^ = J2x'^°'{^)- this definition C,'^^{i3,L) is an effective correlation length which 
converges to the second moment correlation length ^2 in the limit L 00. In our 
computations we used the particular definition (as e.g. in ref. 0): 



2sin(7r/L)\ Gi^h) 



^2(0) , , 

' 1, (4.9) 



where = (2tt/L, 0). 

In our analysis of the Monte Carlo data we shall make the working assumption that one 
is allowed to replace the limiting procedure lim^^^^ limi^oo by 

g^ = limgniz), z := L/C'' {P, L) , 

gn{z) := lim «,(/?, L). (4.10) 

p—>Pc, z iixed 

That is we attempt to first take the continuum limit at fixed physical volume and after- 
wards take the physical volume to infinity. The 2; ^ 00 limit of ^r(-2) is expected to be 
reached exponentially; for example in the leading order 1/n expansion 

9k{z) = ^r(oo)(i - cy/zexp{-z) + ....) . (4.11) 



The situation may however be slightly more complicated due to our particular definition 
of C''- Indeed in the continuum limit at fixed physical volume we expect G'2(0)/G'2(^o) ~^ 
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G{Q)/G{k) where k ~ Kq = (27rMi5/z,0) and the continuum expressions are in finite 
physical volume. On the other hand for the continuum two point function defined in 
infinite volume 



1 



G(0)/G(i^o) - 1 



Ml 



/2vr\2 / 

(t) (^^ - ' 



(4.12) 



In our simulations the values of 27r/z are ~ 1 i.e. not so small; nevertheless at such values 
the correction factor on the rhs of ( [4.12| ) only deviates from 1 by the order 10~^. This 
deviation is much smaller than the statistical accuracy of our simulations, and hence we 
ignore these additional effects in our analyses of the lattice data. 



5. The Ising model 



The particular field theory we are considering in this section is that obtained from the 
Ising model in zero external fieldQfor < T-T^ ^ 0. The spin-spin correlation functions 
in the scaling limit are known exactly from the work of Wu et al |T^, and from this 



knowledge Sato, Miwa and Jimbo found that the S-matrix operator was given by 

S= (-l)N(N-l)/2^ (5.1) 

where N is the particle number operator. An energy independent phase is not observable 
in a scattering experiment; the non-trivial S-matrix ( |5.1|) refiects the fact that the off-shell 
spin-spin correlation functions are not that of a free field. The continuum limit of the 
Ising model is also described by a free Majorana field, but this is non-local with respect to 
the spin field; for a more detailed discussion we refer the reader to the lectures of McCoy 



*One can obtain an infinite number of field theories from the Ising model in the presence of an external 
field H by taking the hmit H ^0,T ->Tc with h = H/\T - Tc\^^/^ fixed. 
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5.1 Form factor determination 



The generalized form factors are given by 



V 



(2,)(A^-i)/2 Yl Ti\e,-e,\) n r(ft - ^ ^ n (5.2) 

l<i<j<m l<r<m<s<N ^ V*^^ "s) rn<k<l<N 

with an odd (positive) integer. We evaluate the dominant contribution to the coupling 
using Eq. (|H^). The non-integrated part ( |3.31| ) vanishes. For the integral ( p.34| ) we need 
fi{0), which is readily obtained from ( |5.2|) , 



M9) = -2iTi29)/T\e). 
Thus the dominant contribution to 74 is 



(5.3) 



74;121 



27r Jo 



du 



16 



5 
'2 
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(5.4) 



Numerically this gives 74:121 = —4.993427441(1) or ^ 14.98 in the leading approxima- 
tion. 



The simplicity of the form factors ( p.2| ) also makes the Ising model a good testing ground 
for the computation of the sub-leading contributions, to which we turn now. The evalua- 
tion of the spectral moments ( p.39|) , ( |2.43 ) is straightforward. For m = 3, 5, 7 the results 
are given in Table ^ 



m 


72;m 


S2 


m 


3 


8.1446256566(1) x 10"^ 


1.094(1 


) X 10-5 


5 


7.96(1) X 10-^ 


2.22(1) 


X 10-10 


7 


7.8(1) X 10-1° 


4.6(1) 


X 10-15 



Table 1: m-particle contributions to 72, S2 in the Ising model 
Table ^ suggests that the series ( |2.43| ) converge extremely rapidly and we would estimate 

72 = 1 + 8.15259(1) X 10-\ 52 = 1 + 1.094(1) X 10-5 , (5.5) 

where the estimated errors come both from the numerical integration and from estimating 
the contributions of the higher particle terms. To get some check on this we may consider 
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the ratio ^2/72 for which from the leading terms ( p.5D we get ^2/72 = 0.999196 336(11). 
This is in excellent agreement with the result ^2/72 = 0.999196 33 of Campostrini et 
al. ^3], which they obtained by numerical evaluation of the exact formula for the 2-point 
function^ of Wu et al. |l9|. 

The evaluation of 74 is more involved. In order to gain insight into the rate of decay of 
the higher particle contributions as well as their sign pattern we pushed the computation 
up to + / + m < 8. The k + I + m = 8 contributions in particular turned out to 
be a formidable computation despite the deceptive simplicity of the form factors. The 
computation is based on the formulae ( p.9| ), ( |3.56| ) and similar ones for (m, 2, m), with m 
odd, and for (1, /, 3), with / even. To give the reader a chance to follow the computations 
we have collected some intermediate results in Appendices B, C. The final results for the 
contributions of the k-l-m intermediate states with A; + / + m < 8 to 74 are summarized 
in Table ^ 



fc, /, m 


74;fc«m 


1,2,1 


-4.993427441(1) 


1,2,3 


0.046310(1) 


1,4,1 


-0.002653(1) 


3,2,3 


0.0002884(3) 


1,4,3 


-0.0000420(5) 


1,2,5 


0.00002562(2) 


1,6,1 


-0.0000040(1) 



Table 2: A;-/-m-particles contributions to 74 in the Ising model 



The rapid decay of the terms is manifest. Increasing A; + / + m by 2 gives a contribution 
roughly two orders of magnitude smaller than the previous one. The sign pattern appears 
to follow the rule: Sign(74;fcim) = Sign(/c + m — / — 1). Further terms with larger differences 
\k — — m\ are suppressed as compared to those with smaller ones. In view of Table |^ 
we would thus (conservatively) estimate the A; + Z + m > 10 particle contributions to be 
< 10% of the sum of the A; + / + m = 8 contributions. This gives 

74 = -4.90321(3) . (5.6) 

'''this famous Fredhohii determinant (solving the Painleve III equation) is basically the summed up 
FF series; see e.g. 
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Inserting into ( |2.42| ) with 72, 82 taken from ( |5.5| ) then yields our final result 



= 14.6975(1) . (5.7) 

This amounts to a determination of g^^ to within < 0.001%. For comparison we collected 
the results of some previous determinations in Table |^ below. 



Method 


Value for g^^ 


High temperature 


14.6943(17) m, 14.67(5) |2|] 


Borel summation 


15.5(8) H 


Monte Carlo 


14.3(1.0) 



Table 3: Previous determinations of (^r in the Ising model 



Finally we would like to mention that an analogous 4-point coupling can be defined 
at criticality T = Tc by sending the magnetic field H to zero. Of course in this case 
the definition of Binder's cumulant has to be modified appropriately to take into account 
the fact that the field has non-vanishing vacuum expectation value. Remarkably can 
be computed exactly by taking advantage of the fact that the small H behavior of the 
partition function is known exactly [^. The final result is = — = —478.307. 



5.2 Recent Monte Carlo simulation of the Ising model 

Our Monte Carlo investigation of g^^ was performed on several IBM RISC 6000 worksta- 
tions at the Werner-Heisenberg-Institut. 

In this subsection is denoted simply by ^. We studied the dependence on the lattice 
spacing by running at /3 = .418 (^ = 10.839936), (3 = .4276 (^ = 18.924790) and (3 = 
.433345 (^ = 33.873923) on lattices of size L = 80, L = 140 and L = 250, respectively. 
These values were chosen in such a way that they have almost exactly the same value 
of z = L/C, ~ 7.4. Figure |l] shows that there is no significant dependence on the lattice 
spacing (i.e. ^). Therefore we decided to use all the data together to study the finite size 
effects. 

We studied the finite size dependence by measuring in addition g^^ on lattices of size 
L = 40, 60, 80, 140 at (3 = .418, (^ = 10.839936). Finite size scaling works very well, i.e. 
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Figure 1: The data for (yf^ at 2; ~ 7.4 for different lattice spacings 



the results only depend on z = L/^. The dependence on z is still quite well described by 
Eq. ( [4.111) . This can be seen in figure |^. A least square fit produces 



The fit quality is not fantastic (x^ = 2.4 per d.o.f.) but acceptable. So our final Monte 
Carlo estimate for g^^ is 



We report our numbers in Table ^ In this table we also indicate the number of measure- 
ments. These were performed using the cluster algorithm as follows: one run consisted 
of 100,000 clusters used for thermalization, followed by 20,000 sweeps of the lattice used 
for measurements. Each run was repeated after changing the initial configuration. One 
such run was considered as one independent measurement. The error was computed out 
of this sample by using the jack-knife method. 

Our estimated value for g^^ in Eq. ( p.9| ) is in very good agreement with the values from 
the analysis of the high temperature expansion given in Table ^ it is also consistent with 
the value Eq. ( |5.7| ) obtained from the form factor construction. 



c = 3.91(3) 



^r(oo) = 14.69(2) . 



(5.8) 



g^ = 14.69(2) . 



(5.9) 
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Figure 2: The data for g^^ vs y/zexp{—z) together with the fit 



6. The XY-model 



In this section we compute the leading contribution to the four-point couphng in the 
two-dimensional 0(2) nonlinear a-model better known as the XY-model. Starting from 
the lattice formulation, after a chain of mappings consisting of several steps the model 
is transformed to a system equivalent to the two-dimensional Coulomb gas. The contin- 
uum limit of the Coulomb gas model (corresponding to the Kosterlitz-Thouless critical 
point [^]) is thought to have a dual description in terms of a Sine-Gordon model at the 
(extremal) Sine-Gordon coupling = Svr. For a review of the XY-model, see . 

In the following we will start by discussing the XY-model S-matrix. The next step is to 
solve the Smirnov equations for the three-particle form factors, which enter the formula for 
the leading term. A general method for finding the Sine-Gordon form factors is given in 
|]3T| . This extends the results of Smirnov where the form factors for an even number of 
particles were found. The spin three-particle form factor we are interested in is probably 
similar to the three-particle form factor of the fermion operator (corresponding to the 
equivalent massive Thirring- model description), explicitly given in Here however we 
need the three-particle form factor only for special rapidities and we found it simpler to 
obtain this special version by going back to the functional equations. It is then used to 
numerically evaluate the leading contribution to (7^. 
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^ of runs 


teflf 


A 


Or 


. t: X (J 


4n 


200 


1 83QQ36 


1 63 54ri 3'l 


1 1 941 (^ 1 

X X . X I X X J 


41 R 


60 


900 


1 8'^QQ'^6 


1 79 81 ri 1 


1 4 1 04l'96'l 




80 


321 


1 839036 


1 73 94f6'l 


14 587f30') 


.418 


140 


300 


10.839936 


174.08(4) 


14.743(39) 


.4276 


140 


100 


18.924790 


455.34(35) 


14.610(54) 


.4276 


250 


225 


18.924790 


455.90(12) 


14.796(46) 


.433345 


250 


202 


33.873923 


1254.48(72) 


14.567(37) 



Table 4: Ising data for x = 62(0) and g^ {p, L) 



6.1 The XY-model S-matrix 

We will regard the XY-model as the n = 2 member of the family of 0{n) a-models. 
Recall that the formulae ( p.2|) , ( p.4|) have a smooth n ^ 2 limit; this has been noted and 
commented on previously by Woo [Q. In this limit 

^i(^) = 7-^^2(^), a2(^) = 0, a,{e) = S2{e) (6.1) 

and 

S2(^) = -exp{22 — sin(M^2(^)}, K^{u^) = 1^ ■ (6-2) 

In this paper we will assume that the spectrum of the XY-model in the (massive) con- 
tinuum limit consists of an 0(2) doublet of massive particles whose S-matrix is given 
by ( |6.1|) with (|6.2|) . Of course, taking the formal n — 2 limit of the bootstrap results 
valid for n > 3 would not be convincing in itself, but ( |6.1| , |6.2|) actually coincide with 
the 87r limit of the Sine-Gordon S-matrix, the prediction of the Kosterlitz-Thouless 

theory! The S-matrix ( |6.1| ) and the corresponding scattering states as a consequence have 
a Wq=_i(su(2)) Hopf algebra symmetry, which as a Lie algebra is isomorphic to su(2). 
The latter is an explicit symmetry in the alternative chiral Gross-Neveu formulation of 
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6.2 The three particle form factor 

Next we calculate the three-particle form factor at the special rapidities necessary to 
compute the leading contribution (|3.34|) . For this purpose we note that the equations for 
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the functions k, I given in subsection 3.2 can relatively easily be solved in this particular 
case n = 2. We first note that Eq. ( |3.54| ) simplifies 

K{m-e) = K{m + e), 

K{in + 0) = -^{{irr -e){in -2e)L{tn -e) - {tn + e){tn + 2e)L{i7T + 6)} .{6.3) 
Inserting ( p.55| ) yields 

r^/ ■ /IX 311^ — 29 in — 29 . ^, ,^ ^, 

K{-m -9)= ■ —— K{-m + 9) . 6.4 

3z7r + 29 m + 29 

Luckily a term proportional to Kiin — 9) drops out here and one is left with the simple 
form ( |6.4| ). This can easily be converted into the form ( |D.1D and solved as 

K{in -9) = {29 - 5ni){29 - 7ni) e^(i) (pUch-) . (6.5) 

Here 

D{9) = Ai{9) + A3_{9) (6.6) 

in the notation of Appendix D and (j){z) is a polynomial function to be determined later. 

Since Y{9) already has the right singularity structure the functions K{9) and L{9) are 
analytic in the physical strip. The residue axioms determine their value at 6' = and 

^ = f as 

K{0) = 0, ^(0) = 1, 

= ;?(fy %) = -»=(t)%)- 

So far we have established that the solution can be expressed in terms of 

Y{9) = -2^^P- e2A(-+e)+A(2^))-A(0) /g_g) 

sh| ch6' 

and 

K{9) = {29 + 3ni){29 + bni) e^^"^) 0(sh-) . (6.9) 



2^ 

The polynomial (j){z) can be determined using the residue constraints { p.7\) , which we can 
rewrite as 

K{0) = 0, K'{0) = -, 

k{^-^) = e^(o)-2^(f ) , ^( - y) = -e^(°)-^(^)-^(^) . (6.10) 
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Using ( |6.8| ) and ( p.4| ) one sees that for real 6 oo 

\Y{e)\r^e~'ei. (6.11; 



This can be used to infer that the polynomial (j){z) can be at most second order, otherwise 
the integral contribution to the leading term would diverge. Taking into account that 
K{0) = and the requirement of real analyticity one must have 

=i0l2 + 02^', (6.12) 



for real constants 0i and 02- Now it is easy to see that ( |6.7| ) determines 0i as 

^ e-^(^). (6.13) 



157r3 



In order to determine 02 we employ the following identities ||33[| 



00 



II chfavro;) — 1 ^,,1 an 1 
pia) := exp<^ / duo \ , \ 6""'^^ = 7 , 6.14 



, f°° , ch(a7rci;) — 1 ,^,,1 1 — , , 

ga := exp / dc^ \. \ e-^"" = ^ , 6.15 

^ sh(7r^) J cos(f) 







to obtain 



exp { A(0) - 2A(|) . - } . Ml) g ^ ^ . (6.16) 

This can be used to show that (|6.10|) is satisfied for the choice (j)2 = 0. Thus </'(sh|) = 
#ish|, and since 461(9) = (in - 2e)K{e) - (in + 2e)K{-e), both k{e) = Y{e)K{e) and 
l{6) = Y{6)L{6) are known explicitly for the XY model. 



6.3 Calculation of the leading contribution 

Having all the ingredients at our disposal we can compute the leading term ( ^.71 ) of the 
intrinsic coupling. Firstly from ( p.32| ) we have for n = 2 

71521 = -(In 4-1). (6.17) 
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Further substituting the exphcit results for the functions A;, / obtained above into Eq. ( |3.41| ) 
and evaluating the resulting expression numerically we obtain 



7i3i = -5.14902(1) , (6.18) 

and hence 

74;i2i = 7g2i + l^lli = -4.65718 . (6.19) 
Thus the leading contribution to the XY-model four-point coupling is 

= -2^ ^ -274;i2i = 9.314 . (6.20) 

I2O2 

Since '^2^2 > 1 and since the next leading contributions to 74 are probably positive (as 
they are in the Ising and 0(3) models), we expect that the true value of g^^ will be less 
than that given in ( 6.20|) (probably by 2 — 4%). 



6.4 Comparison with lattice results 

For the XY model with standard action Kim p| gives the value 

= 8.89(20) (6.21) 

for (3 = 1/0.98. We are in the process of producing higher precision Monte Carlo data 
for this model; so far we can only give a preliminary result, obtained on a lattice of size 
L = 500 at /5 = 1.0174: 

(7k = 9.14(12), (6.22) 

We will return to this issue in a separate publication, where we intend to analyze the 
finite size corrections as well as the lattice artifacts. 

We also wish to mention the results from the high temperature expansion: Butera and 
Comi fT^ obtain 

(7k = 9.15(10), (6.23) 
whereas Pelissetto and Vicari give 

(7k = 9.01(5). (6.24) 

So there is an overall rough agreement between the lattice and the form factor results, 
but the precision is not comparable to that obtained for the Ising model. 
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7. The 0(3) nonlinear sigma- model 



The 0(3) nonlinear sigma model is an important testing ground for quantum field theo- 
retical scenarios in nonabelian gauge theories. The form factor technique has been par- 
ticularly fruitful in studying its possible off-shell dynamics and can be confronted with 



what can be achieved by perturbation theory or numerical simulations |^4[. The intrinsic 
coupling has been computed before by a number of different techniques; we compare the 
results with ours at the end of this section. The present form factor determination takes 
as usual the Zamolodchikov two-particle S-matrix |^ as its starting point; it is given by 



Eqs. (|0|, ^ with n = 3 and 

9 — 711 



S2{e) = -—^. (7.1) 



711 



The corresponding kernel ( |3.4| ) is simply given by 

k^{u) = e-^^. (7.2) 



7.1 Form factor determination of g^i 

Following the by now routine procedure we first collect the ingredients for the evaluation 
of the dominant (1,2, 1) contribution to the intrinsic coupling. From ( p.32|) one readily 
finds for n = 3 

71521 = - • (7.3) 

TT 

The 0(3) form factors have been computed in |^, ^ In particular the reduced 

3-particle form factor Q in Eq. (|3.21| ) is given by 

(7.4) 

where 

rAr(01,...,M = n <d,-e,), 
l<i<j<N 

T^{(^ — ivr) 
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Correspondingly the functions k, I parametrizing via Eq. ( p.40[ ) are for n = 3 explicitly 
given by 

29 7c^T\29)9{29 -7ii) 



m 



m 



m-u AT\9){ti^ + 9^f 
Plugging this into the general formula Eq. ( p. 411) yields 



ill) _ 1 p , r7r6M2(4M2 + 7r2)(2M2 + 7r2) T^{2u) 64 
74;i2i-— /_ dn| ^ _____ 



(7.6) 



57r Jo 



T^{u)ch.^u 



(7.7) 



Numerically we then obtain 74;i2i = —4.16835492(1), so that as a first approximation 
^ — |74;i2i = 6.9472. This is already in rough agreement with other determinations 
in the continuum theory: the 1/n, the e- and the (/-expansions [jl3|, 0. The leading order 
1/n computations have been performed in [36|. For the spectral integrals the result is 



11 11 
72 = 1 + 0.00671941- + 0f—) , ^2 = 1 + 0.00026836- + Of—) . (7.8) 

n ^n^' n ^n^^ 

and for the coupling [|10| 

Stt r 1 / 1 s 1 

(7.9) 



87r 

9r = — 
n 



0.602033- + 0f^ 

n ^n^ 



which gives the approximation g^^ ^6.70 for the case n = 3. The results from the other 
methods are given in Table ^ Considering the rather short series in each case it is amazing 
how well the estimates by the various methods agree. 

For a more precise determination we now return to the form factor approach and examine 



the sub-leading contributions. Using the exact form factors ||3J] the results for the 3- and 
5-particle contributions to 72 and 62 are readily evaluated and are listed in Table ^ 



m 


72;m 


^2;™ 


3 


1.67995(1) X 10-=^ 


3.46494(1) X 10-^ 


5 


6.622(1) X 10-^ 


7.114(1) X 10"^ 



Table 5: m-particle contribution to 72,(^2 in the 0(3) model 

The size of the higher particle contributions to 72 and 62 can roughly be estimated by an 
off hand extrapolation of Table ||; essentially they are negligible to the desired accuracy. 
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The latter could also be justified by referring to a more refined extrapolation scheme, 
based on the scaling hypothesis of ref. In upshot we obtain 



72 = 1.00 1 687(1), 52 = 1.000 034 657(1). (7.10) 

The computation of the sub-leading terms to 74 is much more involved. The starting 
point is again the formulae ( p.lO|) in subsection 3.1. Due to the complexity of the form 



factors however the computation is feasible only computer aided. The essential steps are 
given in appendices B,C. The computation has been performed independently by subsets 
of the authors using slightly different techniques. The final results for the contributions 
of the k-l-m intermediate states with + / + m < 6 to 74 are listed in Table 



A;, Z, m 




1,2,1 


-4.16835492(1) 


1,2,3 


0.051748(1) 


1,4,1 


-0.004065(1) 



Table 6: fc-Z-m-particles contribution to 74 in the 0(3) model 

The leading 1-2-1 contribution is a factor ~ 42 greater in magnitude than the sum of k-l- 
m contributions with k + I + m = 6. It is difficult to bound the rest of the contributions, 
especially since the signs appear to be alternating. The computation of the states with 
I + m + n = 8 would be quite an undertaking. But assuming that the pattern in Table |^ 
continues, as it seems to be the case in the Ising model (see Table §), then we consider 
the assumption that the sum of the remaining contributions fc + Z + m>8is< 10% of 
the sum of the k + I + m = 6 contributions to be reasonable and we then obtain 

74 = -4.069(10) , (7.11) 

and hence our final result 

(7, = 6.770(17) . (7.12) 

This amounts to a determination of to within 0.3%. For comparison we give some 
results of other already published determinations in Table |^. The first two are continuum 
methods while the last one is based on the lattice regularization. We describe the two 
lattice techniques in somewhat more detail in the next subsection, including in particular 
our own recent Monte Carlo results. 
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Method 


Value for g^^ 


5f-expansion 
e-expansion 
High temperature 


6.66(6) § 
6.55(8) 11 
6.56(4) irg, 6.6(1) [11] 



Table 7: Other determinations of gR in the 0(3) model 



7.2 Lattice computations of 

High temperature expansion: 

The analyses of the high temperature expansion for the spectral moments give 72 = 
1.0013(2) |3§ and 62 = 1.000029(5) [0. The agreement with the FFB values Eq. (|7l0|) 



is acceptable; note that these are smaller than that anticipated from the leading order of 
the 1/n approximation, Eqs. ([7.8|). 

The various Fade approximations show the coupling falling rapidly as (3 increases in the 
region of small f3, then a region of rather flat behavior after which these approximations 
show diverse behavior; some analyses indicate that in fact there is a shallow minimum 
and that the continuum limit is actually approached from below (see e.g. refs. |r^). 
In ref. [jll| Campostrini et al. quote for the case n = 3 the result g^ = 6.6(1), and in a 



more recent publication Felissetto and Vicari cite 6.56(4) Butera and Comi on the 



other hand are rather cautious, and did not quote a value for the case n = 3 in ref. |T2 
if pressed they would at present cite g^^ = 6.6(2) ^8 . 

Numerical simulations: 

Monte Carlo computations of g^^ have a long history, see e.g. refs. J^, |^. In order to 
attempt to match the apparent precision attained in the FFB approach, we recently 
performed new high-precision measurements. These were performed on several IBM RISC 
6000 workstations at the Werner- Heisenberg-Institut. In addition we made use of the 
SGI 2000 machine of the University of Arizona, especially for the very time consuming 
simulations on large lattices. 



Based on the fixed point action |T8[ we have measured g^^ at three different values of f5: 
0.70, 0.85 and 1.00, corresponding to correlation length ^ ^ 3.2, 6.0 and 12.2, at the 
values of z = L/^ in the range 5.4 ... 8.2. The data and their analysis can be found in 0, 
the final result is g^^ = 6.77(2). 

Monte Carlo measurements with the standard action were performed using a method 
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367 


34.500(15) 


1267.20(57) 


6.665(14) 


6.730(14) 


1.8 


500 


382 


64.790(26) 


3838.76(1.50) 


6.691(15) 


6.733(15) 


1.9 


910 


127 


122.330(74) 


11883.0(6.4) 


6.737(21) 


6.792(21) 


1.95 


1230 


68 


167.71(17) 


20901.4(19.0) 


6.792(40) 


6.853(40) 



Table 8: 0(3) data for ^, x = 36*2(0) and 5, 



similar to the cluster estimator of |^ . We have reported the analysis of such simulations 
already in our earlier paper . But in the meantime we produced more data and we take 
the opportunity to report them here. 

The present status of the results of our simulations are given in Table In this table we 
also indicate the number of measurements. These were performed using the cluster algo- 
rithm as follows: one run consisted of 100,000 clusters used for thermalization, followed by 
20,000 sweeps of the lattice used for measurements. Each run was repeated after changing 
the initial configuration. One such run was considered as one independent measurement. 
The error was computed out of this sample by using the jack-knife method. 

Our measurements were taken at 6 different correlation lengths ranging from about 11 to 
about 168 on lattices satisfying L/C, ^ 7. To study the finite volume effects, we took in 
addition data at ^ ~ 11 for lattices of sizes L with L/C ^ 5.5, 9 and 13. As discussed 
in the finite size effects are well described by the formula ( |4.11 ), even at finite (large) 
correlation lengths. In the 0(3) model the n = 00 value c = v^Svr fits very well. 

But unlike the Ising model, the lattice artifacts are by no means negligible. To study 
them, we first use Eq. ( [1.1 1| ) to extrapolate our data to z = 00. In this extrapolation we 
use the effective correlation length ,^eff and neglect the fact that this is not exactly equal 
to the exponential correlation length. In Fig. ^ we plot those extrapolated values of g^i 
against 1/^ which we identify with l/C^^. 

Unfortunately there is no rigorous result concerning the nature of the approach to the con- 
tinuum limit. At the time of our last analysis the data point at the largest value of ^ ~ 
168 was not available. In that paper we fitted the data in the entire range from ^ ~ 11 to 
C ~ 122 using a Symanzik type ansatz of the form 5'r(0 = 5'r.(oo) [1 + ^1^^^ log'C + ^2'C~^], 
and thereby obtained the result (^r = 6.77(2). When we now repeat the same fit for the 
new data, which in particular includes the new point at ^ ~ 168, the result is only slightly 
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Figure 3: The extrapolated values of QrI^, oo) 

changed to g^^ = 6.78(2) but the quality of the fit becomes poorer. The fact that the two 
data points closest to the continuum limit lie above 6.78 is in this scenario interpreted as 
a statistical fluctuation. 

On the other hand the present rather large central value at ,^ ~ 168 could be interpreted 
as an indication that the continuum limit is approached much slower than conventionally 
assumed, perhaps as slow as 1/ln^ (which may be expected in the 0(2) model [0)! If 
we adopt this viewpoint it is clear that, although qualitative fits can be made, without 
further analytic information, our data are not sufficient to make a reliable quantitative 
extrapolation to the continuum limit. However, independent of the assumed form of the 
approach to the continuum limit, if the large value at ^ ~ 168 is confirmed by more 
extensive studies it would practically establish a discrepancy between the form factor 
and the lattice constructions of the 0(3) sigma-model. This point, which needs complete 
control over all systematic effects, albeit extremely difficult on such large lattices, is 
certainly worthy of further investigations. 



8. Conclusions 

A new technique to compute the intrinsic 4-point coupling in a large class of two-dimensional 



43 



QFTs has been developed and tested. Starting from the form factor resolution of the 
4-point function the termwise zero momentum limit turned out to exist, providing a de- 
composition of the coupling into terms with a definite number [k, I, m) of intermediate 
particles. Based on the exactly known form factors these terms can be computed prac- 
tically exactly and in the models mainly considered (Ising and 0(3)) were found to be 
rapidly decaying with increasing particle numbers. There is every reason to expect that 
this trend continues, which allowed us to equip the results with an intrinsic error estimate. 
The final results are 

Ising model: (?r = 14.6975(1) , 
0(3) model: = 6.770(17) . (8.1) 

They amount to a determination of g^^ to within < 0.001% and 0.3%, respectively. 

In addition we obtained the universal, model-independent formula ( |A.3| ) for the dominant 



contribution to the coupling, which typically seems to account for about 98% of the full 
answer. We illustrated its use in testing our proposed bootstrap description of the XY- 
model. It would surely also be interesting to apply it e.g. to supersymmetric theories, 
where alternative techniques are hardly available. 

The comparison with the lattice determinations of g^i is quite impressive in the case of 
the Ising model, where there is also very good agreement between the high temperature 
and the Monte Carlo determinations. For 0(2) we are so far lacking both precise Monte 
Carlo and form factor data, but at this preliminary stage there is rough agreement. We 
intend to return to this model in a separate publication. 

The situation in 0(3) is not completely clear: There is a less than perfect agreement 
between the high temperature result and the new high precision Monte Carlo data, and 
there is also room for doubt about the agreement between Monte Carlo and form factor. 
We cannot resolve this question at the moment, mainly because even with our enormous 
amount of Monte Carlo data it is at the moment not clear what the correct extrapolation 
to the continuum is. 
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A. General formula for the dominant term 



Here we describe the generalization of the formula ( p.7[ ) for the dominant 1-2-1 particle 
contribution to g^i to general integrable QFTs without bound states and operators other 
than the 'fundamental' field. The latter is particularly natural in the form factor approach 
because 'fundamental' and 'composite' operators are treated on an equal footing. Thus 
let Oi be possibly distinct, possibly non-scalar but parity odd operators Oi and write oi 



for the quantum numbers labeling them. Parallel to ( 2.13 ) we define the Green functions 
by 

...,kL) = (27r)2<5(2)(A;i + . . . + k^) G'°i-°^(A;i, . . . , A;^) , (A.l) 

where S°^'"^^{ki, . . . , ki) is the Fourier transform of the connected part of the Euclidean 
correlation function {Oi{xi) . . . Ol{xl))- The obvious generalization of the intrinsic cou- 
pling is 

(^01020304 1'n n n fTi 

(A.2) 



E,<fcG°^°'^- (0,0)2 



Here M is again the mass gap and the constant M is conveniently adjusted to nor- 
malize the 1-particle contribution to the denominator to unity. If JF° are the constant 
1-particle form factors of the leading 1-particle contribution to G°^°'^{Q,Q) is just 



^o-i02 ._ ]Vf~2jF°iC"^jF^^, where is the charge conjugation matrix associated with the 
given S-matrix (c.f. below). Thus we take M = '}2j<k{'^°^°^Y ■ With these normalizations 
the dominant 1-2-1 particle contribution to the coupling (|A.2| ) is 



1-2-1 
1 



^ i-oo Aoi ^ 

Jo An 



+ 



16ch^U ' " ' " ~ ~ ' 030201 \- ; -7 -/- - - 030201 \ ' -/ J • (^-3) 

Here the symmetrization is over all elements of the permutation group ^4. D is defined 
in terms of the given bootstrap S-matrix S'^{6) by 



D 



oi02;o304 



'T'Ol 'T'O'2 /^cc' /^dd -^03 -^04 



(A.4) 



and 6*2, 6*3) is the 3-particle form factor of O. Taking the results for the 0(n) 

models as a guideline one would expect that ( |A.3| ) typically yields about 98% of the full 
answer for the coupling. 



In the following we describe the derivation of ( |A.3| ). In contrast to that of ( p.7| ) we keep 
track here of the distributional terms like (|2.19|) and show explicitly that they cancel out 
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in the final answer. In particular this illustrates that the use of the simplifying limit 
procedure ( p.23| ) is justified. 



To fix conventions we first recall the defining relations of a generic bootstrap S-matrix. 
A matrix- valued meromorphic function S^^{6), 6' G (D, is called a two particle 5'-matrix if 
it satisfies the following set of equations. First the Yang-Baxter equation 



STiOi2)S'/M3)St,i023) = Sri023)SZi0i3)S'^i{9^,) , (A.5) 

where 6*12 = 61—62 etc. Second unitarity (|A.6| a,b) and crossing invariance (|A.6| c) 

ST'iO) S:L{-0) = Si^t (A.6a) 

SZ{6) S^^{2m -6)= 5t5t (A.6b) 

Stt{6) = Caa'C"' Sliiir^ - 6) , (A.6c) 



where ( |A.6| c) together with one of the unitarity conditions ( |A.6| a), ( |A.6| b) implies the 
other. Further real analyticity and Bose symmetry 

Finally the normalization condition 

s'am = -s:5t . (A.8) 

The indices a,b, . . . refer to a basis in a finite dimensional vector space V. Indices can 
be raised and lowered by means of the (constant, symmetric, positive definite) 'charge 
conjugation matrix' Cab and its inverse satisfying CadC^^ = ^a- The S-matrix is a 
meromorphic function of 6. Bound state poles, if any, are situated on the imaginary axis in 
the so-called physical strip < Im^ < tt. From crossing invariance and the normalization 
(|A.8| ) one infers that S^^^in) = —CabC"^ is always regular, in contrast to S'^{—m) which 



may be singular. 

Next we prepare the counterparts of Eqs ( p. 101) , ( |2.11 



{^mk){k\0,\l^{l\0,\m){m\OM « Qr'"'{^m. (A.9) 
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where 

niT'^^^m = c^''^''j^T,,{e\o) . (a.io) 

From the S-matrix exchange relations it follows that I^°^{6) is a completely symmetric 
function in 6^ = (6*1, . . . ,9^)- Similarly Iklm"^"^ symmetric in each of the sets 
of variables u = (tui, . . . , Uk), C, = . . . and 9 = {9i, . . . , 9^)- As before we denote 
by Vm{ki,k2) and Vkim{ki, k2, k^, k^) the quantities ( p.8| ) with the integrations over the 
rapidities performed, where the measure is inherited from (|2.25|) . For simplicity we drop 
the operator labels Oj in the notation. When evaluated at kj = {Msh.Kj,0) we write 

Vklm{K.l, K,3, ^,4) , etC. 

In the next step one inserts the expressions for the generalized form factors in terms of 
the ordinary form factors; see for an account in the present conventions. For m = 2 
one obtains explicitly 

.3, = 8ch^.,ch^.4 ich.,(ch., + ch.,)^('^^ + '^^^ 



ch/t2(chKi + chK2) 



Jo Att 4ch^| + (shKi + shKs)^ 
^^Zn{ - /«i + «7r - ze, A, - ^, A, + ^)^°|;(fi:4 + ivr - ze, A, - |, A, + |)*| . 

The Ki —>■ limit of this expression can be evaluated on general grounds. The key 
observation is that a three particle form factor has the following universal 'small rapidity' 
expansion 

1 1 



J'"iO, + t7r-te,92,9s) 



6*12 - ie 6*13 - ie 



d 



+2CaA0isD:^i, - OuDZ,)^°4 , where D^t = -^^S:ti9) 



(A.12) 



9=0 

This expression is uniquely determined by the following properties: (i) The numerator 
is linear and boost invariant in the rapidities, (ii) It obeys the (linearized) S-matrix 
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exchange relations in 62 and 9^. (iii) It has simple poles at ^21 + *e and + ie with 
residues dictated by the form factor 'residue equation' (see e.g. |T^ for an account in the 
present conventions). Using ([A.12|) in (|A.11|) one can compute the small Ki behavior of 



■^121- We denote by V121 the contribution from the non-integrated part and by V121 that 
from the integrated part. One finds 



V^^^ = ^Z"^°-'Z"-'"^[5{Ki + Kg) - 5(Ki + K^)] - 1Z}"1°2'°304 



(A.13) 



where '=' again indicates that both sides give the same result for the symmetrized 
limit and L)°i°2'°304 = £)02ovmo3 = 1)0^03,020^ jg gj^g^^ j^y 

Analyzing the small Ki behavior of the integral in (|A.11|) (by splitting it according to 
du = Jq du + du, e 0^) one finds a distributional term and a regular one. The 
result is 



r°° du 
An 



1 



16ch u 



jro.j,o,^aa3^tb3j,02^^^^ (m, -U, u)^^^'^ C^^'^ J^^^,^,^ ( 



^^0104 ^0203 _|_ ^0103 ^020i^ 



(A.14) 



where the integrand is regular for -u — >■ 0. 



For the generalization of the Q term ( |2.19| ) one obtains in the 1-particle approximation 
and in the k, — ^ limit 



'{h, k2, h, h] 



kj = {MshKj,0) 



5{ki + K2)Z°'"'Z 



03 04 



(A.15) 



Finally, combining (|A.13|) , (|A.14|) and ([A.15|) according to ( |2.21| ) one sees that, - as 
promised in section 2.2 - all distributional terms drop out when computing the right 
hand side of ( |2.21| ). The final result thus does not depend on any prescription how to 
take the Ki ^ limit and is given by ( |A.3| ), as asserted. 
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B. Computation of the 1—4 — 1 contribution 



We start from the general formula Eq. (|3.9| ) with / = 4. For (3i > (32 > P^ > P^ we have 
for one of the factors occurring in ( p.lO| ) 

(a, a\S\<d)\h,(3i- 62, (32] h, (3^- 64, (3^^ = 
Kb^b2b;bS^^ + a - /32, /33, (3i) 

+ATl[5abA0^ - Pl)J'b2b,b,W2,P3,PA) 

+5{a - P3)^debSf^l^ P2, (3i)Sb2bz,ef{(32 - P3)Sbif,daiPl " Ps) 

{P2-PA)Sb,h,da{Pl-PA)]. (B.l) 

For the 5-particle form factor we introduce the reduced form factor through 



( 



h, 02,03, 04,05) = T^{0\, ■ ■ ■ ,05)Qaia2a3a4a5iOl, 02,03, 04,05) 



(B.2) 



Multiplying out we obtain 



'141 



(B.3) 



where we momentarily omit the arguments {k^, —Ki, Pi, P2, P3, Pa)- The shorthands are: 



^{i) ^ ^ (Fl^^b2b3bS^T^ -til- ie, Pi,p2, p3. Pa) 

X ^ib^b2b;bS™ + - ie, Pi,P2, p3. Pa)* , 
)IC^''\k4,-ki,Pi,P2,Ps) 

+SiP4 - K4)IC^''\~Ki,K4,Pl,P2,P3)*} 

HPa ^ P3) + {Pa ^ P2,P2 ^ P3,P3 ^ Pa) 
+ {P4 ^ Pi, Pi ^ P2, P2 ^ P3, P3 ^ Pa) , 

IC^^^n = ^4^f)ciin)^^^^^^^^p^^p^)^s{Pi + Ki)6{P2-K4) + {Pi 
HP2 ^ P3) + {P2 - Pa, Pa ^ P3, P3 ^ P2) 
+ {Pi ^ P2, P2 ^ P3, P3 ^ Pi) + iPi ^ P2, P2 ^ Pa, Pa 
+ (/3i ^ /33, P3 ^ /32, /32 ^ Pa, Pa ^ A) , 



/?3,/33 



(B.4) 



(B.5) 



A) 



(B.6) 
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where 

K:^''\a, 7, A, /32, /33) C •^M.bslA, /32, /33):^iV6.63(^^ + a - ie, 7, A, /32, /^s)* (B.7) 



61,62,63 



IC^^^^\k4, -Ki, (33, (5^) =^ .^616364 (/«4, /?3, /54)^626364(-«;i, /?3, /?4)*'5l62,6il(«;i + «;4)(.B.8) 

6162,63,64 

Because of the symmetry in the I3i arguments of 1C^^\1C^^^\1C^^^^^ one has 

, J.. ]^ pco poo poo poo 

iZZOoTT J— 00 J— 00 J— 00 J— 00 

X —4 — '{K4,-Ki,Pi,P2,P3,P4) , (B.9) 

"J poo poo p 00 ^ 

^^^l?('«l>'«2,/^3,/^4) = =7^-^ / / d/^a / d/jg „ — (B.IO) 

7687r"' J-oc J-00 J -00 1 + Efc=i ch/Sfe 

.Kf)(«:„«:2,«:3,«:4) = ^ f d/?i f d/32^^^^^^^£(^^^)(«^4, A, /32) • (B.U) 

D47r 7 -00 z + 2^k=i ciipfe 

The contribution (///) is very simple; we can set the Ki to zero to obtain 
where we decomposed 

(0,0, /?,-/?) = T^{2(5)T\(5)k^"'\(5) . (B.13) 

Writing similarly 

K:^''\a, 7, A, /32, /33) = k^''\o^, 7, A, /^s, Pz) 



- 7) i<^<j<3 k=i T{l3k -Oi-ie) 

one has 

^i4?=i:^£'^\ (B.15) 
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with 



(II 1] 11 f°° 

(«1,«2,«3,«4) = l^^^T^f \ V / ^/^l / / d/33 

oo J —oo J —oo 



X- 



7687r2 T{ki + K4) 
n nPi-Pj) 



1 + Efe=i ch/3fc i<f<-^.<3 



fc=l 



/?2, /?3, -«:3)fc^''^ /^4, A, /92, /Ss)* n T(/?fc - t^,)-—^——} , (B.16) 

fc=l ^ [Pk + Kl) ^ 



Vul K2, K3, K4) = -— / d(3i / d/52 

izoTT J-00 J-00 



k=l 



1287r J-00 J-00 2 + Efe=i ch/?fc 

xr2(/3i - /?2) n T{Pk - i^A)T{(5k + ACi) 
fc=l 

X {fc(^^)(«^4, -/^i, /3i, /?2, «:4) - «;4, A, /32, • (B.17) 

In the latter term we can set the Ki to zero to obtain 

(B.18) 

Lastly we turn to the (J) contribution. There arc many ways to manipulate the integral 
into a form more amenable to numerical evaluation. Here we proceed as follows: Writing 

/C(^) = A;«K-«i,/3i,/32,/33,/?4) n T'{Pi-Pj) 

l<j<j<4 

^ 1 

H r(/5, + + te)T{(3k -k,- te) ' ^^"^^^ 

we replace the l/(x ± ie) distributions by a sum of products of principal parts and delta 
functions, thereby obtaining 

3 



The three terms are: 

Vl^l'iKi, K2, K3, K4) 



122887r3 



d/ii / d/52 / d/33 / d/34- 



C» J — c» 
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A. A, ft. A) ^ J^^T=(ft - ft) n . (B.21) 



15367r2 + K4) 

/oo foo poo 1 

-00 ^-00 J-00 i- + 2^k=l^'^h>k i<i<j<3 

X {6{Pu /?2, /33, «:2)A;(^)(«:4, A, P,, P,, -k,) {[ T(/?fe + k,)—-^ 

+S{(3,, (32, (3s, -m)k'-'\K,, -K,,(3,, (32, (3s, k,) {[ T{(3u - ^4) ^.^^ J (B.22) 



C7 3") "J- /""^ 

V\^i'{ki, K2,K3, K4) = — — / dPi / d/3; 

ZODTT J-00 J-00 



''2 + ELich/?fe 
2 

k^'\K4, -Ki, A, /32, K4)T\P, - P2) n - «4)T(/3fe + «i) . (B.23) 



k=l 



In the latter expression we can set the Ki to zero to obtain 

«llf'(o, 0.0,0) = IZ •i^ d.pii + .^.^f ^mrWHo, 0, ft -ft 0. 0) . (b.24) 

For W^(^) := 

■^^li'i^^ we now invoke the identity 

ni<.<,<4 sh(y. - y) ^ 1 f (-l)^ch(y,) sh(,.-,,), (B.25) 

to get 

W^^^ = Jini^ [W^^^ [A] {ac) + W^^^ [B] (aj] , (B.26) 

with the notation 

W^'^[X]{ac) = 12 daiGxiai) , X^A,B. (B.27) 

Here 

'P roo 'p roo 

Ga(q;i) = T\a) J '^^'^ T{a ) J '^'^3-^a(q;i, ^2, as) 
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Gs(q;i) = 

where 

In these formulae 
Tx{o(i, a2, as) = 



y roc y roo ^ 

1 {ai) Jo I [a2) Jo 

-^a(q;i, -a2, as) - ^a(-Q;i, 02, as) + ^^(-ai, -0^2, as) •} , 

G{ai) _ 4G(0) 
^h^ 

Q/-. roo POO 

G(ai) = ch^— / da2 da3J^ij(ai, a2, "s) • 

Z J —00 J —00 



1 

16 



/c(^)(0, 0, «i, q;2, "3, «4)/x(«l, "2, "3, «4) IlLi ^h^^ 



cha4(E^=i cha^) Yli<j ch^^^ 

where 7 is given through sh7 = — J2k=i shafc and 

/A(Q;i,a2,a3,a4) = -3sh sh 

X sh sh sn~ 

2 2 2 

, 2Q;2 - as , 2Q;2 - "4 , 2Q;3 - «4 

jB{ai,a2,a3,a4) = sh sh sh . 

For the [B] contribution it is numerically convenient to decompose 

W^^\B]{ac) - hG{0) + W^^^BO] + W'^^^[Bl]{ac) , 



04=7 



where 



and 



h 



Finally we recombine 



sh^f ' 

- dai \^ 4t1 I = -0.0014539754 . 



3 

Vui{Ki, K2, K3, K4) = W^^\ki, K2, K3, K4) , 
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with 



w(i) _ ,,(^.1) 



W^'^ = v\^^^+v^'\ (B.38) 



Case n = 1: 

Here we simply have (recall the 2-particle S-matrix = — 1) 

fc(^) = 16 , = -8i , = -4 , (B.39) 

from which one sees 

^ilf (0,0,0,0) = -t;il{'')(0,0,0,0), 

(0,0,0,0) = v[[{'\0,0,0,0) = -^v[i{''\0,0,0,0) . (B.40) 

Thus 

W^^'^ = = W^^'^ , (B.41) 

so that for the Ising case we simply get Vi^i = W^^\ with W^-^^ given by Eq. ( [B.21| ) and 
by Eq. (p39| ). This is, as expected, the same expression as that obtained directly 
with the form factor written as a product over principal parts as in Eq. ( p.2[ ). 

Case n = 3: 

Firstly for W^'^^ we obtain 

= -0.0005420(1) . (B.42) 

Next for W^"^^ one has 

1 1 roo roo roo 

W^'\t,,,t,2,^3,^d = ^ , / d/3i / d/?2 / d/53J(/?l,/32,/33) 

00 J— 00 J—OD 

^ fc=i T{/3k - K4) 

+S{(3,,(32,f33, ^3)w^'\^^i, ^A,(3,,(32,f33) n + ^^^)7f7^ , (B.43) 

k=i T(/3fc-Ki)J 

where 

J(/?i,/?2,/33) = — n (B.44) 
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and 

k^''\K^, -Ki, A, /?2, /?3) + lk^'\l^4, -1^1, A, /?2, /?3, -/^l) , (B.45) 
W^'^\ki, Pi, P2, Ps) = 

-A;(^^)(-«:i, K^, -f3i, -f32, -^3)* + '^k^'\fii, -^1, -f3i, -f32, -/^a, Hi) • (B.46) 
Explicit calculation reveals the fact 

^^^'H«i,«4,/5i,/32,/53) =«;^'H«4,«i,/3i,/32,/53). (B.47) 
Now we expand w^^^ for small 

W^'^\ki, Ki, Pl, p2, Ps) = WoiPi, p2, Ps) 

+ {ki + K4)Wi{pi, p2, Ps) + {ki - K4)W2{P1, p2, Ps) + ^(k^) . (B.48) 

In fact we do not require W2- Note that the functions Wi are real, so that in particular 

lmw^^\0,0,Pi,P2,P3)^0, (B.49) 
which is needed to avoid a singularity in W^'^^ for ^ 0. Hence 

= (2) [A] + W^(') [B] , (B.50) 

with 

1 1"°° r'^ 1 

^^^^ " 192^ /-oo'^'^' y_^d/32^J(/?i,/32,/3o)«'i(/3i,/32,/3o) , (B.51) 

where /3o is determined through sh/3o = — sh/3i — sh/32 , and 

11 POO fOO fOO 

x{5(/3i,/32,/33,«;2) n^(/^fe + '^i)™ 7 + {1^2 ^ Ki ^ 1^4)} , (B.52) 

fe=l T{Pk-Ki) > 



with 

^(A, /i^2, P3) = J{Pi, P2, P3)wo{Pi, P2, P3) ■ (B.53) 
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We then see that Vr*^^^[-B] is a sum of two parts 



with 



I fOO roo 1 "T? 

,,..IB2] ^ -1^ /" d,. dft J^A (MM) 



Numerically this gives 



W^'^\A] = 4.41085(1) X 10"^ 
W^^^Bl] = 4.9600(1) X 10^^ 
W^^\B2] = 1.1503(1) X 10-^ 



and hence 



= 0.00050219(1) . 

Finally we turn to the computation of W^^\ Due to Eq. (p.49| ) it follows that 

t;r4i''^(0,0,0,0) = -2i;ilf(0,0,0,0). 

Now explicit computation yields 

k^'\0,0,P,-P,0,0) = 7r-6|r3(0,/?,-/5)P(40/?2 + 327r2), 

So for W^^^ we arrive at 



64 7o ch^/3 /36(/32 + 47r2)^ 

-0.0004682756 . 
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C. Computation of the 1 — 2 — 3 contribution 



We use the results (|3 .571 - ^1611) with m = 3, and begin with contribution {IV): 



1 /■ 3 6{ai,a2,a3,-K4) 



15367r2 J Y: .=i chai 

X Q''^^\-K2 + ^vr_, -Ki + i7r_, «!, ^2, as) • (C.l) 

Here 

where A stands for 9i, 62, 63, 64^, 65 and A' for 6^3, 6*4, 6*5. Note that 5f('^^)(y4) is totally 
symmetric in the subset A' . 

We decompose the 1/ (x ± ie) factors to obtain terms involving products of principle parts 
and delta-functions; only terms having less than three delta-functions contribute in the 
Lim procedure i.e. 



■'123'' — ^'^123' ^ ■ (C-3) 



s=l 

The terms are 

(/y,i). N -L . . f ,3 S{ai,a2, 03,-1^4) 



VuVHf^i, 1^2, K3, 1^4) ~ i53g7j-2 ^(^i ~ '^2) J d^a- 



xc/(^^)(m, OT, ai, aa, "3) H ^^("^ " "i) 11 ^ V^F? \ V ' 

^ /"j2 5( — /«2, ttl, ^2, — '«4) 



Vl23'^^('«l,'«2,K3,K4) ~ --7^ / d^a 



1 + Ej=i chai 



X5(('^^)(i7r,i7r,0, ai,a2)T^(ai - 0L2)\\T[ak ^ 1^2)— ■ -, (C.5) 

while for the term involving the product of two delta-functions one finds 

^m'^H'^i? ^2, «3, '^4) = 0{k^) . (C.6) 

The contribution (/V, 1) is antisymmetric in k\ ^ K2 and so it doesn't contribute in the 
sum over permutations. In the contribution (/V", 2) we can take the k — limit and 
obtain 



57 



Now turn to the (///) contribution: 



where 

G^^"\Oi, 62, 03,64) = X] •^16162(^1' ^2, 6'3)^ft\b2i(^2,6'3, 614)* 

6162 

= 73(^1, 02, ^3)^3(^2, ^3, o,rg^'''\e,, 02, e,, e,) . (c.9) 

Here one can set the Hi to zero to obtain 

^{23'^ (0, 0, 0, 0) — r d.l3^X^ -.9^'''^ (^TT, 13, -13, 0) . (C.IO) 

123 V , , , y 12877 7-00 ^ch^/3[l + 2ch/3]^ v y \ ) 

For the (//) contribution: 

(//)/ X 1 f .2jiPl,P2,Ki,K2) 

J ch[3i + 2^1=1 chctj 

where 

= T3(A')T3(5)T3(5')*^^''n^) , (C.12) 

where A stands for ^1, ^2, ^3, ^4, ^5; ^' stands for ^1, ^2, ^3, -B stands for ^3 + i7r_, ^4, ^5, and 
stands for 62,04,95. Making the famihar decomposition of the singular distributions 
one obtains: 

4? = E43'^\ (C.13) 

s=l 



with 



,(11,1) 



I \ 1 /'^2^ W^/^2,^1,^2) if^... 

i^u n., ns, -^)-^2j^ P ehA + ch/32 ^^-^'^ 



/ d^ce '^f^^'"^'^^' /^'^^ 5f^^^^(-/ti +CT,/3i,/32,Qii,a2)J'^(Qii - ^2) 

- n TTvT^ n m -«.■); ^ 
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V,2, i^U ^2, .3, ^4) - y d /3 ^^^^^ ^ ^^^^j ^ ^^^^ ^ ^^^^j (C.15) 

+ ZTT, /32, (32, -Ks)T\(3, - ^2) H T{l3j + /^a) ^ 



f^S'^Vi, /^s, «;3, ^4) ~ y" d^aT^(ai - ^2) (C.16) 



V 



——2 r T '[-Ki+n:,-K2,-Ki,ai,a2)- 

1 + Ei=i chtti 



In the contributions s = 2, 3 we can set the Ki to zero to obtain 

«^'(o.o,o,o) ~ ^/;;d^-pi!M_,<")(„,^,-A-^,o), (C.17) 

«*(o, 0,0.0) ~ i^/;;d. ^j;^^;^^j ,'"'(...o.o.„,-„). (ds) 

Finally for the (/) contribution 



+ ch/32 



/ d^a^i^^|ai^l^^W(-«i + i7r_, A, /32, Ci, ^2, a,) , (C.19) 



with 



6 01,02,03 

= r3(>i')7^5(s)T3(s')V'H^)- (C.21) 

Here A stands for 61, 62, 6*3, 64^, 9^, 9q; A' stands for 9i, 92, 6*3, B stands for 9^ + 6*2 + 
m, 9i,9^,9Q, and B' stands for 9^,9^,9^. 

Then we rearrange to terms where after doing the (52 integral the singularities in the /3i 
integral all have negative imaginary parts 

4^3 = E4^- (C.22) 

s=l 



59 



The terms are: 

(7,1). V 1 f 2j{Pu(^2,Ki,K2) 

/ d a —3 ^(-Ki + z7r,/^i,/i^2,Q;i,Q;2,a3) 

T{Ki + Pi + i€)T{Ki + P2-i€)IlkT{Pi + ie-ak)T{P2-i€-ak)' ^ ' ^ 



ch/?i + chp2 



/ ^ ^ 1 ;o , v^2 T—r 'y-i^i + «7r, /32, ai, q;2) 

^ m-/^2)T^K-«2)n.m-^») 24) 

+ a + ie)r(fi:i + /?2 - ie) Hi r(/?2 - ie - a^-) ' 

^'m Vi> '«2, /t3, «4) ~ 30 727^2 ^*-^^ ~ '^2'' 

/ d a —3 g"^ ^(-Ki + i-K, -K2, -Ki, ai, a2, aa) 

X— n^<^-^'K;-%) (C.25) 

rifc T{-Ki -le- ak)T{-K2 + le - ak) 

.g?K,«:2,.3,«:4)~-^TK-.2) 

f -,2 S{-K2,ai,a2,-K4} . 

/ d a — —^^2 r 9 '[-i^i + ^tt, -^2, -«2, «!, a;2) 

.^ ^^("l -«2)^^^(ft:2 + «^) .^26) 

As for the last contribution, it vanishes in the hmit Ki —>■ 

4^(0,0,0,0) = 0. (C.27) 



For the contribution (/, 1) we now perform the P2 integral and shift the /3i integral to 
larger imaginary part, after which we can send all the to zero to obtain 

1 r+oo+i<l> _ch^^ 



«-°-°) -512^ -^V/? 

roo roo 

/ dui du2T\ui)T\u2)T\ui+U2)M^^\uy 
Jo Jo 
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xg(^\zn, (3, ai, a,, as) > (^-28) 

where the ak are determined in terms of the m's as in Eq. (p.30|) . 
For the (/, 2) term we obtain 

viif -vgf+vgf + 0{K), (C.29) 

where 



.,3 (-1, -2, -4) - ^ y d ,hA + ,h/32 

d a— ——2 7—^^ ^(-Ki + ITT, pi, Pi, ai, a2) 

chpi + Ei=i cha^ 

xT(ft - A)T=(a. - a.) R ^j;^ 11 ^(ft - -^)^^^) • (C.30) 

.,23 (-1, -2, ^s, n,) --j^jd /^ ^^ehA[l + E,ch/j.] 

+ ^vr, A, /?2, A, -P2, -^^3)T\Pi - P2) n T{P^ + -(C.Sl) 
In the latter we can do the (32 integral and set the Hi to zero to obtain 

(0. 0. 0. 0) ~ £ ^■^ oh-.a^fwl ^''''"- ^' °' ■ 

Finally 

=S'(0. 0. 0. 0) = - jl^ £ d„-j-J(M_,<',,,,, 0, 0, „, 0) . (C.33) 

Summarizing the results we have 

1^123(0, 0,0,0) = ^\/(^), (C.34) 
i=i 

where the five terms cirG clS follows: 

y«=43'^(0, 0,0,0), (C.35) 
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given in Eq. ( |C2(j| ). Further 



1 .nT^{W) 



d/^^l^^^'H/^) , (C.36) 



where 

(2),., _ 9^'''\^^^, -P, 0) - tg^'^'KtTr, tTT, 0, (3, -/3)ch/? 

Next 

= t;!2?(0,0,0,0) + t;f//)(0, 0,0,0) 

where 



(7(3) ai, aa) = 2(7(^^)(z7r, /5, -/?, ai, ^2) + ^^^^'^m, /5, /5, ai, ^2) . (C.39) 



Further 



where 



Finally 



where 



V'"' = i.i'2f(0,0,0,0)+iS"(0, 0,0,0) 

1 ,c.4o) 



2567r j-00 cha[l + 2cha]' 



^(4)(«) = -,^(3)(o,a,-a). (C.41) 



0^) = t;f2f(0,0,0,0) + t;S'^(0, 0,0,0) 
1 r T\2(3) 



/(5)(/?) = -2(7(3)(/?,-/?,0). (C.43) 
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Case n — 1: 

Here we have simply 

g'-'^A) = -16, g^''\B) = 8i, 



and so 



Thus 



(7(')(/3)=4, ^W = 0,r = 3,4,5. 
sh^/? 1 



with 



327ry-oo ch^/3 487r ' 

- -T2^f <'^"^"(-)^^(-)^^("'+-)il0■ 



Numerically we find 

= -0.000842721(1) 

Case n = 3: 

First we have 



^^'^ = 32768i, d'«lV^M^(«2)V'K + «2)rM(3n«)-''5(«i,«2), 



9 , 9 



where 
and 

S{ui,U2)^ / d/3 -g- — /i^ ^(iTT, /?,-/?, q;i,q;2, as) 

J-oo+i4> ch p 

(/32 + ^2)3 11 _ /32)2 + 27r2(a| + /32) + Vchafc - ch/3^ 
Numerically this gives 

1^(1) = -0.000844527(1) . 
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Doing the contractions yields 

3 -BO) = -_l(ll±l!ll!m 
g^'''\t7r,t7T,0,(3,-P) = 2tg^'''\tn,P,-(3,0). (C.53) 

Thus 

g'^'\P)=g'^'''\zn,P,-P,0), (C.54) 

and 

512 7-00 '^/52ch^/3(/32 + 7r2)2(/52 + 47r2) 
Next by exphcit computation one verifies 



T/(2) r sh'/^ (2/32 + 7r2)(4/j2 + vr2) 

^ - 512 loo ^^^M^^T^^TWf^^Ti;^ = 0.0074380765 . (C.55) 



g^'\a) = 0, gi^)((3) = -gi^)(-(3)^ (C.56) 

and thusQ 

= v^{5) = . (C.57) 
It remains to compute V^^\ Shifting the ai integral we obtain the representation 
= df3 da- ^ ' 



2567r2i_oo "^i-oo shf ch=^/3 T(a) 
cha,ch^[ch/3 + cha, + cha/ ^"^ " "^^^^"^ " ^^^^"^ " ^^^^T '^^-^'^ 

where 

«!=«-/?, (C.59) 

and a2 is determined by 

shQ;2 = — sh/5 — shai . (C.60) 

Numerically this gives 

= -0.000125112(1) . (C.61) 



Eq. ( p. 57 ) is perhaps true for all n but we have not verified this conjecture 
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D. Building blocks of form factors 



In form factor calculations one often encounters the problem of finding an analytic function 
f{9) satisfying 

f{e) = aie)fi-e), 

f{m-e) = f{in + 9), (D.l) 
for given (7(9). If a (9) has the Fourier representation 

rOO r\(.) ^ 

ai9) = e'^(') , 6(6) = 2 — sm{9ij) k{uj) , (D.2) 

JO UJ 

with some kernel function k{ijj) then the 'minimal' solution of ( |D.1| ) is given by 

The function A(^) has the following properties. If k{(jj) ~ e~^^ [z > 0) for a; — > oo then 
A(^) is analytic for —z < lm6 < 27r + z and for real 6 ^ oo 

Re A{e) ~ A(m + e)r^-- k{0) - — k'{0) + const. (D.4) 

2 vr 

We encountered in Section 6 the following special case: for some (positive, real) parameter 
a 

Mg)-e-M"'^ ;;^°'" + ^ . (D.5) 

(1 + a)in — 

corresponding to the kernel 

Kico) = -e-""(^+") . (D.6) 
We denote the corresponding solution by Aq,(^^). 
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